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Abstract. This is a survey of recent results on nodal sets of eigenfunctions of the Laplacian 
on Ricmannian manifolds. The emphasis is on complex nodal sets of analytic continuations 
of eigenfunctions. 
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Let {M,g) be a (usually compact) Riemannian manifold of dimension n, and let {(fj} 
denote an ortlionormal basis of eigenfunctions of its Laplacian, 

(1) Ag ifj = -A| v^j {^j, Vk) = 5jk. 

Here {u,v) = j^^uvdVg where dVg is the volume form of {M,g). If dM 7^ we impose 
Dirichlet or Neumann boundary conditions. When {M,g) is compact, the spectrum of A 
is discrete and can be put in non-decreasing order Aq < Ai < A2 t oo- The eigenvalues A^ 
are often termed energies while their square roots Xj are often termed the frequencies. The 
nodal set of an eigenfunction cpx is the zero set 

(2) Z^;Q={xeM:Mx) = 0}. 

The aim of this survey is to review some recent results on the "H"^ ^-surface measure and on 
the yet more difficult problem of the spatial distribution of the nodal sets, i.e. the behavior 
of the integrals 

(3) f / fdSx^, {feC{M)) 

as A — cxD. Here, dSx = rfH""^ denotes the Riemannian hypersurface volume form on 
Z^p^. More generally, we consider the same problems for any level set 

(4) ■■= = c}, 

where c is a constant (which in general may depend on A). Nodal sets are special level sets 
and much more attention has been devoted to them than other level sets, but it is often of 
interest to study general level sets and in particular 'high level' sets or excursion sets. 

Date: March 3, 2013. 

Research partially supported by NSF grant # DMS-0904252. 

^In difference references we use either the notation Z oi J\f for the nodal set. Sometimes we use the 
subscript Lp\ and sometimes only A. 
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We have recently written surveys |Z5tlZ6j on the global harmonic analysis of eigenfunctions, 
which include some discussion of nodal sets and critical point sets. To the extent possible, 
we hope to avoid repeating what is written there, but inevitably there will be some overlap. 
We refer there and jH] for background on well-established results. We also decided to cover 
some results of research in progress (especially from [Z3j, but also on quantum ergodic 
theory). We generally refer to the results as 'Conjectures' even when detailed arguments 
exist, since they have not yet been carefully examined by others. 

There are two basic intuitions underlying many of the conjectures and results on eigen- 
functions: 

• Eigenfunctions of Ag-eigenvalue — are similar to polynomials of degree A. In 
particular, Z\ is similar to a real algebraic variety of degree A. 

Of course, this intuition is most reliable when (M, g) is real analytic. It is quite 
unclear at this time how reliable it is for general C°° metrics, although there are some 
recent improvements on volumes and equidistribution in the smooth case. 

• High frequency behavior of eigenfunctions reflects the dynamics of the geodesic flow 
G* : S*M S*M of M. Here, S*M is the unit co-sphere bundle of {M,g). 

When the dynamics is "chaotic" (highly ergodic), then eigenfunctions are de- 
localized and behave like Gaussian random waves of almost fixed frequency. This 
motivates the study of Gaussian random wave models for eigenfunctions, and suggests 
that in the 'chaotic case' nodal sets should be asympotically uniformly distributed. 

When G* is completely integrable, model eigenfunctions are highly localized and 
their nodal sets are often exhibit quite regular patterns. The latter heuristic is not 
necessarily expected when there exist high multiplicities, as for rational flat tori, and 
then some weaker randomness can enter. 

Both of these general intuitions lead to predictions about nodal sets and critical point sets. 
Most of the predictions are well beyond current or forseeable techniques to settle. A principal 
theme of this survey is that the analogues of such 'wild' predictions can sometimes be proved 
for real analytic (M, g) if one analytically continues eigenfunctions to the complexification 
of M and studies complex nodal sets instead of real ones. 

As with algebraic varieties, nodal sets in the real analytic case are better behaved in the 
complex domain than the real domain. That is, zero sets of analytic continuations of eigen- 
functions to the complexification of M behave like complex algebraic varieties and also reflect 
the dynamics of the geodesic flow. 

It is well-known that the complexification of M can be identified with a neighborhood of 
the zero-section of the phase space T*M. That is one reason why dynamics of the geodesic 
flow has greater impact on the complex nodal set. 

We will exhibit a number of relatively recent results (some unpublished elsewhere) which 
justify this viewpoint: 



Theorem 8^, which shows that complex methods can be used to give upper bounds 
on the number of nodal components of Dirichlet or Neumann eigenfunctions which 
"touch the boundary" of a real analytic plane domain. 



EIGENFUNCTIONS AND NODAL SETS 



3 



Theorem 9A on the hmit distribution of the normahzed currents of integration 
over the complex zero sets of "ergodic eigenfunctions" in the complex domain. 



Theorem 11.2 and Corollary 11.1, which show that the similar currents for analytic 
continuations of "Riemannian random waves" tend to the same limit almost surely. 
Thus, the prediction that zero sets of ergodic eigenfunctions agrees with that of 
random waves is correct in the complex domain. 

Sharper results on the distribution of intersections points of nodal sets and geodesies 



on complexified real analytic surfaces (Theorem 10.1). 

Our analysis of nodal sets in the complex domain is based on the use of complex Fourier 
integral techniques (i.e. generalized Paley- Wiener theory). The principal tools are the ana- 
lytic continuation of the Poisson-wave kernel and the Szego kernel in the complex domain. 
They become Fourier integral operators with complex phase and with wave fronts along the 
complexified geodesic flow. One can read off the growth properties of complexified eigen- 
functions from mapping properties of such operators. Log moduli of complexified spectral 
projectors are asymptotically extremal pluri-subharmonic functions for all {M,g). These 
ideas are the basis of the articles \Z2\ \TZ\ [Z3t [Zlf [Z8| \Z9\ IHe] . Such ideas have antecedents 
in work of S. Bernstein, Baouendi- Goulaouic, and Donnelly-Fefferman, Guillemin, F.H. Lin 
(among others) . 

We note that the focus on complex nodal sets only makes sense for real analytic {M,g). 
It is possible that one can study "almost analytic extensions" of eigenfunctions for general 
C°° metrics in a similar spirit, but this is just a speculation and certain key methods break 
down when g is not real analytic. Hence the results in the C°° case are much less precise 
than in the real analytic case. 

It should also be mentioned that much work on eigenfunctions concerns ground states, i.e. 
the first and second eigenfunctions. Unfortunately, we do not have the space or expertise to 
review the results on ground states in this survey. For a sample we refer to |Mej . Further, 
many if not all of the techniques and results surveyed here have generalizations to Schrodinger 
operators —h'^A + V. For the sake of brevity we confine the discussion to the Laplacian. 



0.1. Notation. The first notational issue is whether to choose to be the positive or 
negative Laplacian. The traditional choice 




makes is negative, but many authors call — A^ the Laplacian to avoid the minus signs. 
Also, the metric g is often fixed and is dropped from the notation. 

A less traditional choice is to denote eigenvalues by rather than A. It is a common 
convention in microlocal analysis and so we adopt it here. But we warn that A is often used 
to denote A-eigenvalues as is |DF| [H] . 

We sometimes denote eigenfunctions of eigenvalue — A^ by ip\ when we only wish to em- 
phasize the corresponding eigenvalue and do not need ip\ to be part of an orthonormal basis. 
For instance, when A^ has multiplicities as on the standard sphere or rational torus, there are 
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many possible orthonormal bases. But estimates on "H" ^(^i^a) ^'^ depend on whether 
(fx is included in the orthonormal basis. 

0.2. Acknowledgements. Thanks to C. D. Sogge and B. Shiffman for helpful comments 
on the exposition, and to S. Dyatlov for a stimulating discussion of L°° quantum ergodicity. 



1. Basic estimates of eigenfunctions 

We start by collecting some classical elliptic estimates and their applications to eigenfunc- 
tions. 

First, the general Sobolev estimate: Let w G C^{Q) where C with n > 3. Then 
there exists C > 0: 



Next, we recall the Bernstein gradient estimates: 

Theorem 1.1. |DF3j Local eigenfunctions of a Riemannian manifold satisfy: 
(1) Bernstein estimate: 



(6) / |V<^a|W < — / |<^a|W 



B(p,r) / \JB{p,r) 



(2) L°° Bernstein estimate: There exists K > so that 
(7) max |Vv?A(a;)| < max 

xeB(p,r) r xeB(p,r) 

(3) Dong's improved bound: 

max V(j9a < max 09a 

Br{p) r Br{p) 

for T < CaA-i/^ 

Another well-known estimate is the doubling estimate: 

Theorem 1.2. (Donnelly-Fefferman, Lin) and [H] (Lemma 6.1.1) Let ip\ he a global eigen- 
function of a C°° (M, g) there exists C = C{M, g) and tq such that for < r < tq, 

Further, 

c\ 



/ J" \ OA 

max iv^Afa^)! < ( — ) max \Lpxix)\, (0 < r' < r). 

B(p,r) \r / xeB{p,r') 



The doubling estimates imply the vanishing order estimates. Let a E M and suppose that 
u{a) = 0. By the vanishing order z/(m, a) of m at a is meant the largest positive integer such 
that D^uia) = for all la I < z/. 
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Theorem 1.3. Suppose that M is compact and of dimension n. Then there exist constants 
C{n),C2{n) depending only on the dimension such that the the vanishing order v{u^a) of u 
at a E M satisfies z/(m, a) < C(n) iV(0, 1) + C2(n) for all a e i?i/4(0). In the case of a global 
eigenf unction, i'{ipx,a) < C{M,g)\. 

We now recall quantitative lower bound estimates. They follow from doubling estimates 
and also from Carleman inequalities. 

Theorem 1.4. Suppose that M is compact and that ipx is a global eigenfunction, A(^a — 
X^(fix. Then for all p, r, there exist C,C' > so that 

max \(pxix)\ > C'e^'^^. 

xeB(p,r) 

Local lower bounds on |log|</?^| follow from doubling estimates. They imply that there 
exists A,S > so that, for any e Mr/2, 

(9) sup |v^a(C)I >Ce-^^ 

CeB5(Co) 

Indeed, there of course exists a point Xq G M so that |v?A(a;o)| ^ 1- Any point of Mr/2 can 
be linked to this point by a smooth curve ofuniformly bounded length. We then choose S 
sufficiently small so that the 5-tube around the curve lies in M^- and link ^^(C) to Bs{xo) by 
a chain of 5-balls in Mj- where the number of links in the chain is uniformly bounded above 
as C varies in M,-. If the balls are denoted Bj we have sup^^.^^ \(fx\ < e^^sup^^. \(fx\ since 
-Bj+i C 2Bj. The growth estimate implies that for any ball B, sup2B \(fix\ < e'^'^sup^ \ 
Since the number of balls is uniformly bounded, 

1 < sup |<^a| < e'^^ sup \ipx\ 

Bs{xo) BsiO 

and we get a contradiction if no such A exists. 

As an illustration, Gaussian beams such as highest weight spherical harmonics decay at 
a rate e"^'^'^^^'''') away from a stable eUiptic orbit 7. Hence if the closure of an open set is 
disjoint from 7, one has a uniform exponential decay rate which saturate the lower bounds. 

We now recall sup-norm estimates of eigenfunctions which follow from the local Weyl law: 

^x{x, x) Ea.<a WM? = (27r)-" J^^,^^)^, di + i?(A, x) 
with uniform remainder bounds 

|i?(A,x)| ^ CA"-\ xeM. 

Since the integral in the local Weyl law is a continuous function of A and since the spectrum 
of the Laplacian is discrete, this immediately gives 

k.(x)p^2CA-i 

which in turn yields 

(10) ll^^ll^o = 0(A^) 
on any compact Riemannian manifold. 
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1.1. estimates. The classical Sogge estimates state that, for any compact Riemannian 
manifold of dimension n, we have 



11) ^^ = 0(A^(p)), 2<:p^oo, 



where 

(12) 6{p) 



2(n+l) ^ ^ 
n— 1 ^ ^ 
2(n+l) 



2 ^ p ^ 



n-l 



Since we often use surfaces as an illustrantion, we note that in dimension 2 one has for A ^ 1, 



(13) WfxhriM) ^CX'^^-^ p'^yxh^M), 2^p^6, 
and 

(14) \Wx\\lp{m) ^ CX "-^ v' H|eA||L2(M), e^p^cx). 



These estimates are also sharp for the round sphere S . The first estimate, (13), is saturated 



by highest weight spherical harmonics. The second estimate, (14), is sharp due to the zonal 



functions on S'^, which concentrate at points. We go over these examples in ^3.2 



2. Volume and equidistribution problems on nodal sets and level sets 

We begin the survey by stating some of the principal problems an results regarding nodal 
sets and more general level sets. Some of the problems are intentionally stated in vague 
terms that admit a number of rigorous formulations. 

2.1. Hypersurface areas of nodal sets. One of the principal problems on nodal sets is 
to measure their hypersurface volume. In the real analytic case, Donnelly-Fefferman ( jDFj 
(see also jLinj ) ) proved: 

Theorem 2.1. Let {M,g) be a compact real analytic Riemannian manifold, with or without 
boundary. Then there exist Ci,C2 depending only on {M,g) such that 

ciA < W''-\Z^J < C2X, {A^x = AVa; ci, C2>0). 

The bounds were conjectured by S. T. Yau |Yll IY2j for all C°° {M,g), but this remains 
an open problem. The lower bound was proved for all C°° metrics for surfaces, i.e. for n = 2 
by Briining |Br] . For general C°° metrics the sharp upper and lower bounds are not known, 
although there has been some recent progress that we consider below. 

The nodal hypersurface bounds are consistent with the heuristic that (p\ is the analogue 
on a Riemannian manifold of a polynomial of degree A, since the hypersurface volume of a 
real algebraic variety is bounded by its degree. 

2.2. Equidistribution of nodal sets in the real domain. The equidistribution problem 
for nodal sets is to study the behavior of the integrals ^ of general continuous functions / 
over the nodal set. Here, we normalize the delta-function on the nodal set by the conjectured 



surface volume of ^2.1 More precisely: 



Problem Find the weak* limits of the family of measures { j-rf^A }. 
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Note that in the C°° case we do not even know if this family has uniformly bounded 
mass. The high-frequency limit is the semi-classical limit and generally signals increasing 
complexity in the 'topography' of eigenfunctions. 

Heuristics from quantum chaos suggests that eigenfunctions of quantum chaotic systems 
should behave like random waves. The random wave model is defined and studied in 1241 



(see 1:11) , and it is proved (see Theorem 11.1) that if one picks a random sequence {V'Aj } of 
random waves of increasing frequency, then almost surely 

(15) f / fdS^^ ^ — ^ I fdV,, 

A, Jn^^^ Vol{M) Jm 

i.e. their nodal sets become equidistributed with respect to the volume form on M. Hence 
the heuristic principle leads to the conjecture that nodal sets of eigenfunctions of quantum 
chaotic systems should become equidistributed according to the volume form. 

The conjecture for eigenfunctions (rather than random waves) is far beyond any current 
techniques and serves mainly as inspiration for studies of equidistribution of nodal sets. 

A yet more speculative conjecture in quantum chaosis that the nodal sets should tend 
to CLEq curves in critical percolation. CLE refers to conformal loop ensembles, which are 
closed curves related to SLE curves. As above, this problem is motivated by a comparision 



to random waves, but for these the problem is also completely open. In ^|T2 
heuristic principles which started in condensed matter physics |KHt IKHS 



we review the 
isil ITsiKl IWei] 

before migrating to quantum chaos |BS[ IBS2[ IFGSl IBGSj [5^ lEGJSj . It is dubious that such 
speculative conjectures can be studied rigorously in the forseeable future, but we include 
them to expose the reader to the questions that are relevant to physicists. 

2.3. norms and nodal sets. Besides nodal sets it is of much current interest to study 
L'P norms of eigenfunctions globally on (M, g) and also of their restrictions to submanifolds. 
In fact, recent results show that nodal sets and norms are related. For instance, in ^we 
will use the identity 

(16) \\fx\\L^ = ^f \^fx\dS 



relating the norm of (px to a weighted integral over Z^^ to obtain lower bounds on 
7{"-i(Z^J. See 



Obtaining lower bounds on norms of eigenfunctions is closely related to finding upper 
bounds on norms. The current bounds are nowhere near sharp enough to improve nodal 
set bounds. 

2.4. Critical points and values. A closely related problem in the 'topography' of Laplace 
eigenfunctions ipx is to determine the asymptotic distribution of their critical points 

Ciifx) = {x : Vifxix) = 0}. 

This problem is analogous to that of measuring the hypersurface area 'H"'^^{Zx) of the nodal 
(zero) set of (px, but it is yet more complicated due to the instability of the critical point set as 
the metric varies. For a generic metric, all eigenfunctions are Morse functions and the critical 
point set is discrete. One may ask to count the number of critical points asymptotically as 
A — > oo. But there exist metrics (such as the flat metric on the torus, or the round metric 
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on the sphere) for which the eigenfunctions have critical manifolds rather than points. To 
get around this obstruction, we change the problem from counting critical points to counting 
critical values 

CVi^y) = {^^{x) : V^x{x) = 0}. 

Since a real analytic function on a compact real analytic manifold has only finitely many 
critical values, eigenfunctions of real analytic Riemannian manifolds (M, g) have only finitely 
many critical values and we can ask to count them. Moreover for generic real analytic metrics, 
all eigenfunctions are Morse functions and there exists precisely one critical point for each 
critical value. Thus, in the generic situation, counting critical values is equivalent to counting 
critical points. To our knowledge, there are no results on this problem, although it is possible 
to bound the "H"" ^-measure of C{lp\) (see Theorem |Baj ). However 'H"'~^(C((/9a)) = in the 
generic case and in special cases where it is not zero the method is almost identical to bounds 
on the nodal set. Thus, such results bypass all of the difficulties in counting critical values. 
We will present one new (unpublished) result which generalizes eqreflD to critical points. 
But the resulting identity is much more complicated than for zeros. 

Singular points are critical points which occur on the nodal sets. We recall (see [HI IHHLt 
IHHONj ) that the the singular set 

S(^a) = {xe Z^^ : Wipx{x) = 0} 

satisfies 'HJ'~''^{T,{lp\)) < oo. Thus, outside of a codimension one subset, Z^p^ is a smooth 
manifold, and the Riemannian surface measurt 

refer to |HHUNl El [HHLl iHg] for background. 



manifold, and the Riemannian surface measure dS = l v^^ dVg on Z^p^ is well-defined. We 



2.5. Inradius. It is known that in dimension two, the minimal possible area of a nodal 
domain of a Euclidean eigenf unction is 7r('^)^. This follows from the two-dimensional Faber- 
Krahn inequality, 

\k{VL)kma.{D) = Ai(L')Area(L)) => nj^ 

where D is a nodal domain in Q. In higher dimensions, the Faber-Krahn inequality shows 
that on any Riemannian manifold the volume of any nodal domain is > CA~" [EK] . 

Another size measure of a nodal domain is its inradius rx, i.e. the radius of the largest 
ball contained inside the nodal domain. As can be seen from computer graphics (see e.g. 
|HEJj ). there are a variety of 'types' of nodal components. In |Man3] . Mangoubi proves that 

(17) ^ > r-A > 



A - - A5'=(")(logA)2'^-4' 

where k{n) = — 15n/8 + 1/4; note that eigenvalues in jManj are denoted A while here we 
denote them by A^. In dimension 2, it is known (loc.cit.) that 

(18) ^>.>5^. 

2.6. Decompositions of M with respect to (px- There are two natural decompositions 
(partitions) of M associated to an eigenfunction (or any smooth function). 

(i) Nodal domain decomposition. 
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First is the decomposition of M into nodal domains of ipx. As in [PS] we denote the 
collection of nodal domains by ^(v^a) and denote a nodal domain by A. Thus, 

M\Z^,= U A 

When is a regular value of (fx the level sets are smooth hj^ersurfaces and one can ask how 
many components of occur, how many components of the complement, the topological 
types of components or the combinatorics of the set of domains. When is a singular value, 
the nodal set is a singular hypersurface and can be connected but one may ask similar 
questions taking multiplicities of the singular points into account. 
To be precise, let 

/^(v^a) = #^(v2a), i^i'^x) = # components of Z{lpx). 

The best-known problem is to estimate n{fx)- According to the Courant nodal domain 
theorem, /i(v5A,J < n. In the case of spherical harmonics, where many orthonormal bases 
are possible, it is better to estimate the number in terms of the eigenvalue, and the estimate 
has the form ^{<fx) < C{g)\^ where m = dimM and C{g) > is a constant depending on 
g. In dimension 2, Pleijel used the Faber-Krahn theorem to improve the bound to 

limsup — < ^2 < 0.69 

A^oo A Jq 

where jo is the smallest zero of the Jq Bessel function. 

A wide variety of behavior is exhibited by spherical harmonics of degree N . We review the 
definitions below. The even degree harmonics are equivalent to real projective plane curves 
of degree A^. But each point of MP^ corresponds to a pair of points of S'^ and at most one 
component of the nodal set is invariant under the anti-podal map. For other components, 
the anti-podal map takes a component to a disjoint component. Thus there are essentially 
twice the number of components in the nodal set as components of the associated plane 
curve. 

As discussed in [Ley], one has; 

• Harnack's inequality: the number of components of any irreducible real projective 
plane curve is bounded hy g + 1 where g is the genus of the curve. 

• If p is a real projectove plane curve of degree then its genus is given by Noether's 
formula 

_{N-l){N-2) ^ ordp(x)(ordp(s) - 1) 

singular points x 

where ordp(a;) is the order of vanishing of <^a at x. Thus, the number of components 
is < + 1 for a non-singular irreducible plane curve of degree A^. 

Curves which achieve the maximum are called M-curves. Also famous are Harnack curves, 
which are M curves for which there exist three distinct lines of MP^ and three distinct arcs 
aj of the curve on one component so that ^^a^ fl ^j = N. It follows from Pleijel's bound that 
nodal sets of spherical harmonics cannot be maximal for large A^, since half of the Pleijel 
bound is roughly .35A^^ which is below the threshold .5A^^ + 0{N) for maximal curves. 
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Associated to the collection of nodal domains is its incidence graph Fa, which has one 
vertex for each nodal domain, and one edge linking each pair of nodal domains with a 
common boundary component. Here we assume that is a regular value of ip\ so that the 
nodal set is a union of embedded submanifolds. The Euler characteristic of the graph is the 
difference beween the number of nodal domains and nodal components. In the non-singular 
case, one can convert the nodal decomposition into a cell decomposition by attaching a one 
cell between two adjacent components, and then one has fJ^ifx) = '^{fx) + 1 (see Lemma 8 
of IWl)- 

The possible topological types of arrangements of nodal components of spherical harmonics 
is studied in |EJNj . They prove that for any m < N with N — m even and for every set of 
m disjoint closed curves whose union is invariant with respect to the antipodal map, there 
exists an eigenfunction whose nodal set has the topological type of the union of curves. Note 
that these spherical harmonics have relatively few nodal domains compared to the Pleijel 
bound. It is proved in |NSj that random spherical harmonics have aN"^ nodal components 
for some (undetermined) a > 0. 

Morse-Smale decomposition 

For generic metrics, all eigenfunctions are Morse functions [U]. Suppose that / : M — t- M 
is a Morse function. For each critical point p let W^p^ (the stable or descending cell through 
p) denote the union of the gradient flow lines which have p as their initial point, i.e. their 
a-limit point. Then Wp is a cell of dimension Xp = number of negative eigenvalues of Hpf. 
By the Morse-Smale decomposition we mean the decomposition 

M= [j w; 

p:df{p)=0 

It is not a good cell decomposition in general. If we change / to — / we get the decomposition 
into ascending (unstable) cells M = IJp d/(p)=o • intersections fl WgU are always 

transversal then V/ is said to be transversal. In this case dim(PVp fl W^) = Ap — Ai + 1 and 
the number of gradient curves joining two critical points whose Morse index differs by 1 is 
finite. 

We are mainly interested in the stable cells of maximum dimension, i.e. basins of attraction 
of the gradient flow to each local minimum. We then have the partition 

(19) M= U w;. 

p a local min 

This decomposition is somtimes used in condensed matter physics (see e.g. |Wei] ) and 
in computational shape analysis jReu] . In dimension two, the surface is partitioned into 
'polygons' defined by the basins of attraction of the local minima of if. The boundaries of 
these polygons are gradient lines of which emanate from saddle points. The vertices occur 
at local maxima. 
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An eigenfunction is a Neumann eigenfunction in each basin since the boundary is formed 
by integral curves of Vv^a- Possibly it is 'often' the first non-constant Neumann eigenfunction 
(analogously to ipx being the lowest Dirichlet eigenfunction in each nodal domain), but this 
does not seem obvious. Hence it is not clear how to relate the global eigenvalue to the 
Neumann eigenvalues of the basins, which would be useful in understanding the areas or 
diameters of these domains. Note that 

/ <^jdV = / V<^x-udS = 0, 

where z/ is the unit normal to dWp, since V(p\ is tangent to the boundary. In particular, 
the intersection n Wp is non-empty and is a connected hypersurface which separates Wp 
into two components on which (fx has a fixed sign. To our knowledge, there do not exist 
rigorous results bounding the number of local minima from above or below, i.e. there is no 
analogue of the Courant upper bound for the number of local minima basins. It is possible to 
obtain statstical results on the asymptotic expected number of local minima, say for random 
spherical harmonics of degree N. The methods of |DSZj adapt to this problem if one replaces 
holomorphic Szego kernels by spectral projections (see also |Nicj .) Thus, in a statistical sense 
it is much simpler to count the number of "Neumann domains" or Morse-Smale basins than 
to count nodal domains as in |NS] . 

3. Examples 

Before proceeding to rigorous results, we go over a number of explicitly solvable examples. 
Almost by definition, they are highly non-generic and in fact represent the eigenfunctions of 
quantum integrable systems. Aside from being explicitly solvable, the eigenfunctions of this 
section are extremals for a number of problems. 

3.1. Flat tori. The basic real valued eigenfunctions are (pk{x) = sin{k,x) or cos{k,x) 
{k G Z") on the fiat torus T = M^/Z". The zero set consists of the hyperplanes {k,x) = 
mod 27r or in other words (x, |||) G 2^^- Thus the normalized delta function J^dS\z^^ 
tends to uniform distribution along rays in the lattice Z". The lattice arises as the joint 
spectrum of the commuting operators Dj = and is a feature of quantum integrable 
systems. 
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The critical point equation for cos{k, x) is k sin(A;, x) = and is thus the same as the nodal 
equation. In particular, the critical point sets are hypersurfaces in this case. There is just 
one critical value = 1. 

Instead of the square torus we could consider ]R"/L where L C is a lattice of full 
rank. Then the joint spectrum becomes the dual lattice L* and the eigenfunctions are 
cos(/c, x), sin(/c, x) with k G L* . 

The real eigenspace 1-Lx = M — span{sin(fc, x), cos(A;, x) : \k\ = A} is of multiplicity 2 for 
generic L but has unbounded multiplicity in the case of L = Z" and other rational lattices. 
In that case, one may take linear combinations of the basic eigenfunctions and study their 
nodal and critcal point sets. For background, some recent results and further references we 
refer to pZ]. 



3.2. Spherical harmonics on S"^. The spectral decomposition for the Laplacian is the 
orthogonal sum of the spaces of spherical harmonics of degree N, 

oo 

(20) L\S') = ^Vn, A\v^ = XNld. 

N=0 

The eigenvalues are given by A^^ = A^(A^+1) and the multiplicities are given by ttin = 2A^+1. 
A standard basis is given by the (complex valued) spherical harmonics which transform 
by e**"^ under rotations preserving the poles. 

The are complex valued, so we study the nodal sets of their real and imaginary 
parts. They are separable, i.e. factor as CAr^m-P^ (?") sin(m^) (resp. cos(m6') where is 
an associated Legendre function. Thus the nodal sets of these special eigenfunctions form a 
checkerboard pattern that can be explicitly determined from the known behavior of zeros of 
associated Legendre functions. See the first image in the illustration below. 




afsmEIlK)iLM,iK,lKiMtK 



Among the basic spherical harmonics, there are two special ones: the zonal spherical har- 
monics (i.e. the rotationally invariant harmonics) and the highest weight spherical harmon- 
ics. Their nodal sets and intensity plots are graphed in the bottom two images, respectively. 
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Since the zonal spherical harmonics Y^l^ on S"^ are real- valued and rotationally invariant, 
their zero sets consist of a union of circles, i.e. orbits of the rotation action around 

the third axis. It is well known that F^f (r) = ^J^^^Pn 

(cosr), where P/v is the iVth 

Legendre function and the normalizing constant is chosen so that | lyj/^l 12,2(52-) = 1, i.e. 

4:71 J^^'^ \ Pn {cos r)\'^dv{r) = 1, where dv{r) = sinrdr is the polar part of the area form. Thus 
the circles occur at values of r so that P/v(cosr) = 0. All zeros of Pn{x) are real and it 
has N zeros in [—1, 1]. It is classical that the zeros ri, . . . , r^r of P/v (cosr) in (0, vr) become 
uniformly distributed with respect to dr [Si] . It is also known that P/v has — 1 distinct 
critical points [Cl ISz2j and so the critical points of Y^^ is a union of — 1 lattitude circles. 

We now consider real or imaginary parts of highest weight spherical harmonics . Up 
to a scalar multiple, Yn{xi,X2,X2) = {xi + 1x2)^ as a harmonic polynomial on M^. It is 
an example of a Gaussian beams along a closed geodesic 7 (such as exist on equators of 
convex surfaces of revolution). See [R] for background on Gaussian beams on Riemannian 
manifolds. 

The real and imaginary parts are of the form Pj)^(cosr) cos A^6', P^(cosr) sin A^6' where 
P^{x) is a constant multiple of (1— x^)"^/^ so P^(cosr) = (sinr)^. The factors sin A^6', cos NO 
have A^ zeros on (0, 27r). The Legendre funtions satisfy the recursion relation P/^^ = —{2^ + 
\)\/\ — x^P^{x) with Pq = 1 and therefore have no real zeros away from the poles. Thus, 
the nodal set consists of A^ circles of longitude with equally spaced intersections with the 
equator. 

The critical points are solutions of the pair of equations ^P/v (i^) cos N6 = 0, P^ sin A^6' = 
0. Since P^ has no zeros away from the poles, the second equation forces the zeros to occur 
at zeros of sin N9. But then cos N9 7^ so the zeros must occur at the zeros of ^P^{r). The 
critical points only occur when sinr = or cosr = on (0,7r). There are critical points at 
the poles where Y^ vanishes to order A^ and there is a local maximum at the value r = | of 
the equator. Thus, Re has A^ isolated critical points on the equator and multiple critical 
points at the poles. 

We note that |Rey/p is a Gaussian bump with peak along the equator in the radial 
direction. Its radial Gaussian decay implies that it extremely small outside a. N2 tube 
around the equator. The complement of this tube is known in physics as the classically 
forbidden region. We see that the nodal set stretches a long distance into the classically 
forbidden region. This creates problems for nodal estimates since exponentially small values 
(in terms of the eigenvalue) are hard to distinguish from zeros. On the other hand, it has 
only two (highly multiple) critical points away from the equator. 



3.3. Random spherical harmonics and chaotic eigenfunctions. The examples above 
exhibit quite disparate behavior but all are eigenfunctions of quantum integrable systems. 
We do not review the general results in this case but plan to treat this case in an article in 
preparation |Z9j. 

We now contrast the nodal set behavior with that of random spherical harmonics (left) 
and a chaotic billiard domain (the graphics are due to E. J. Heller). 
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4. Lower bounds on hypersurface areas of nodal sets and level sets in 

THE 0°^ CASE 

In this section we review some recent lower bounds on 'H'^~^(Z, 



Before we begin, we recall the co-area formula: Let f : M ^ 
continuous function m on M, 

dV 



from [UMllg^lHSllHW] . 

i be Lipschitz. Then for any 



1 u{x)dV = 




IM J 


M J f 



u- 



df 



dy. 



Equivalently, 



uix)\\Vf\\dV 



M 



ud'H'^-')dy. 



f-Hv) 



We refer to ^ as the "Leray form" on the level set {/ = ?/}. Unlike the Riemannian surface 



measure dS = dTi^ it depends on the choice of defining function /. The surface measures 
are related by dT-L^^^ = |V/|^. For background, see Theorem 1.1 of |HL] . 



Theorem 4.1. As above assume that 



(21) 



A 



2 



1. Then 



M 



for some uniform constant C . Consequently, 



(22) 



A^<H"-^(Z. 



A ^ 1. 



Inequality (22) follows from (21) and the lower bounds in |SoZj 
(23) 



A 4 < 



Iv^aI dVg. 



M 



The lower bound (22) was first proved by Colding and Minicozzi |CMj . A slightly weaker 
result was proved in |SoZj by a different method inspired by the article of R. T. Dong [Dong 
which was successively improved in |HWl IHSj to the same bound as |CMj . 



The L^-lower bounds in (23) were proved in |SoZj using Holder's inequality and the 
eigenfunction estimates of Sogge Sog| for the range where 2 < p ^ ^3^-^^. As will be seen 



n-l 



below, the estimate is sharp for Gaussian beams such as highest weight spherical harmonics. 
The random wave model would predict that Hv^aIIli > C{M,g) > and that would imply 
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Yau's conjectured lower bound in the 'chaotic case'. This provides ample motivation to study 
norms of eigenfunctions on manifolds with ergodic geodesic flow. 

The proof of Theorem 4.1 is based on an identity from |SoZj (inspired by an identity in 
Dong]): 



(24) 



/ \ip^\{^, + \^)fdVg = 2 [ \Vgy^x\f 



dS, 



We recall that dS is the Riemannian surface measure on Z^^ . In |SoZj , we substituted / = 1 
and used (23) to obtain a power lower bound. This was improved in |HWt IHSj by putting 
/ = 1 but applying Schwarz's inequality to get 



/ \^x\dVg^2{'H''-\Z^jY/^ ( [ \Vg^x\^dS 



1/2 



(25) 

The final important step is to show that 



(26) 

by choosing 
(27) 



\Vg^x?dS<X" 



/= (l + AVi + |V,v.A|J)^ 



From (|24|) it follows that 



M 



lifxl (Ag + A^) ( 1 + Xipl + IVgCpxl^) ' dVg. 



By the L^-Sobolev bounds Hv^aHh^^m) = 0{\^) it follows that the "second term" of the right 
side satisfies 



J M 

and thus to prove (26), it suffices to show that the "first terms" satisfies 

V^aI A,( 1 + AVa + \^9^M f'dVg = 0(A3). 



(28 



M 



We refer to [HSj for further details on this bound. A simpler approach was suggested by W. 
Minicozzi, who pointed out that (26) also follows from the identity 



(29) 



\Vqex\^dSq 



M 



sgn(eA) divg( IV^^caI V^ca ) dVg. 



This approach is used in |Arj to generalize the nodal bounds to Dirichlet and Neumann 
eigenfunctions of bounded domains. In the next setction we explain how to obtain more 
general identities. 

There are several ways to prove the identity in (24). One way to see it is that d^x '■ = 
(A + X^)\{px\dV = away from {{px = 0}. Hence this distribution is a positive measure 
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supported on Z^p^. To determine the coefficient of the surface measure dS we calculate the 
limit as 5 —> of the integral 



/(A + A2)|v9,|rf\/ 



M 



/(A + A2)|^aM^. 



IvaI<5 



Here / G C^(M) and with no loss of generality we may assume that 5 is a regular value of 
<y9A (by Sard's theorem). By the Gauss-Green theorem, 

\ /(A + \^)\ipxW - [ Iv^aKA + \')fdV = [ {fd,\^x\ - \y^x\duf)dS. 

Av>x\<S J\'Px\<S A'Px\=S 

Here, u is the outer unit normal and di, is the associated directional derivative. For 6 > 0, 
we have 



(30) 



on {(fx = S}, V 



V<^A 



on {v9a = -5}. 



|v<^aI " |v<^Ai 

Letting 5 — )■ (through the sequence of regular values) we get 

/ /(A + \^)\i^My = lim / /(A + X')\^x\dV = lim / fd^dS. 



Since Iv^aI = ±V^a on {(fx = ±5} and by (30), we see that 



J,J{A + X')\^x\dV = lim,^oL,=,/ 



V\(fx\dS 



1^x1=5 \^\fx\\ 

The Gauss-Green formula and limit are justified by the fact that the singular set S^^^ has 
codimension two. We refer to |SoZ] for further details. 

The lower bound of (23) follows from eigenfunction estimates in Sog , which say that 



I II ^ (n-l)(p-2) 



2<p^?(!Lil). 

n— 1 



If we pick such a 2 < p < ^^^^^p, then by Holder's inequality, we have 

1/0 i-l / (n-l)(p-2) l_i 

which implies Hv^aUli ^ cA~^, since (1 



— (- - -1 



2(P-1)' 



lx (n-l)(p-2) ^ 

4p 4 ■ 



Remark: One can also integrate the identity (24) over a basin of attraction of a local 



minimum (or maximum) (19), since the boundary term vanishes. Thus we get an identity 



between the norm of ipx on each basin and the IVv^aM'S'- measure of the nodal line inside 
the basin. 
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4.1. More general identities. For any function we have 

We then take x to be the meromorphic family of homogeneous distribution x^. We recall 
that for Rea > —1, 

r", X > 



0, X < 0. 



The family extends to a G C as a meromorphic family of distributions with simple poles at 
a = —1, —2, . . . , —k, . . . using the equation ^x+s = sx^~^ to extend it one unit strip at a 
time. One can convert x^ to the holomorphic family 

^+-f(^Ti)' ^+ ■ 

The identity we used above belongs to the family, 
(31) (A + sA2)(^^ = s(s-l)|V^^|Vr'• 



Here = v^*x^ has poles at s = — 1,— 2,---. The calculation in (24) used \(^\ but is 
equivalent to using (31) when s = 1. Then v^+~^ has a pole when s = 1 with residue 



5o{ip) = T^^dS\z ; it is cancelled by the factor s — 1 and we obtain (|24|). This calculation is 



formal because the puUback formulae are only valid when dip ^ when ip = 0, but as above 
they can be justified because the singular set has codimension 2. The right side also has a 
pole at s = and we get Aip^ = —\Vip\'^6'{ip), which is equivalent to the divergence identity 
above. There are further poles at s = —1, —2, . . . but they now occur on both sides of the 
formulae. It is possible that they have further uses. 

The question arises of how such identities are related to the Bernstein-Kashiwara theorem 
that for any real analytic function / one may meromorphically extend to C by construct- 
ing a family Ps{D) of differential operators with analytic coefficients and a meromorphic 
function b{s) so that Ps{D)f^^^ = b{s)f^. In the case / = ipx, the operator | Vv9|~^(A + sA^) 
accomplishes something like this, but it does not have analytic coefficients due to poles at 
the critical points of (f. One wonders what Ps{D),b{s) might be when / = (fx- 

4.2. Other level sets. These results generalize easily to any level set A/"^^ := {(fx = c}. 
Let sgn(x) = ^. 

Proposition 4.2. For any C°° Riemannian manifold, and any f G C{M) we have, 

(32) / f{A + X^)\ipx-c\dV + X^c [ fsgn{(px-c)dV = 2 [ f\V(px\dS. 

Jm J JM^^ 

This identity has similar implications for 'H"~^(A/'^^) and for the equidistribution of level 
sets. Note that if c > sup |v5A(a^)| then indeed both sides are zero. 



Corollary 4.3. For c g 



A^ f ipxdV = I \V^x\dS < X^VoliMf^. 
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Consequently, if c > 

The Corollary follows by integrating A by parts, and by using the identity, 

Im Iv'a - c| + c sgn(^A -c)dV = J^^^^ ipxdV - J^^^^ cpxdV 

(33) 

since = Xv,, ^xdV = J^^^^ ^xdV + J^^^^ ^xdV. 



4.3. Examples. The lower bound of Theorem 4.1 is far from the lower bound conjectured 



by Yau, which by Theorem 2.1| is correct at least in the real analytic case. In this section 



we go over the model examples to understand why the methds are not always getting sharp 
results. 

4.4. Flat tori. We have, | V sin(/c, = cos^(A;, x) |A;p. Since cos(A;, x) = 1 when sin(A;, x) = 
the integral is simply \k\ times the surface volume of the nodal set, which is known to be of 
size \k\. Also, we have | sm{k,x)\dx > C. Thus, our method gives the sharp lower bound 
H"-^(Z^J > CX^ in this example. 

So the upper bound is achieved in this example. Also, we have | sin(/c, > C. 

Thus, our method gives the sharp lower bound ^{""^(Z^^) > CX^ in this example. Since 
cos(/c,x) = 1 when sin(fc,x) = the integral is simply \k\ times the surface volume of the 
nodal set, which is known to be of size \k\. 

4.5. Spherical harmonics on S*^. The L} of norm can be derived from the asymptotics 
of Legendre polynomials 

Pn{cos 6) = V2{7rN sin 9)-^ cos (^{N + - ^) + 0{N-^/^) 

where the remainder is uniform on any interval e < 6 < tt — e. We have 

IICIU^ = ^^V^^^^/^^' \PNicosr)\dvir) ~ Co > 0, 
i.e. the norm is asymptotically a positive constant. Hence IVYq'^M'^ ~ CqA^^. In this 



example |VY'o'^U°° = saturates the sup norm bound. So the estimate of (|23j) produces 



the lower bound "H" ^{Z^J > Aa. The accurate lower bound is A, as one sees from the 
rotational invariance and by the fact that has A^ zeros. The defect in the argument is 
that the bound IVYq^Iloo = A^2 is only obtained on the nodal components near the poles, 
where each component has length ~ 

Gaussian beams 

Gaussian beams are Gaussian shaped lumps which are concentrated on A~^ tubes T^-i (7) 
around closed geodesies and have height A^. We note that their norms decrease 



like A ' 4 ' , i.e. they saturate the L'^ bounds of |Sog| for small p. In such cases we 
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have IVv^'aM'S' ~ A^||(y9A||Li — 4 . It is likely that Gaussian beams are minimiz- 
ers of the norm among L^-normalized eigenfunctions of Riemannian manifolds. Also, 

ri + 1 

the gradient bound ||Vv2a||l°° = 0(A^^) is far off for Gaussian beams, the correct up- 
per bound being X'^"^. If we use these estimates on Hv^aHli and ||Vv9a||l°°, our method 
gives 'H"'~^{Z^^) > CX^~^^, while A is the correct lower bound for Gaussian beams in 
the case of surfaces of revolution (or any real analytic case). The defect is again that 
the gradient estimate is achieved only very close to the closed geodesic of the Gaussian 
beam. Outside of the tube 7^_i (7) of radius A~2 around the geodesic, the Gaussian beam 

and all of its derivatives decay like e~^'^^ where d is the distance to the geodesic. Hence 
IVv^'aM'S' ~ ^j- ^ ^^-j I Vv^'aM'S'. Applying the gradient bound for Gaussian beams to 

the latter integral gives W-^iZ^^ n r 1 (7)) > CA^- V, which is sharp since the intersec- 
tion Z^p^ '"''^-2 (^) '^^^^ across 7 in ~ A equally spaced points (as one sees from the Gaussian 
beam approximation). 

4.6. Non-scarring of nodal sets on {M,g) with ergodic geodesic flow. In this section, 
we prove a rather simple (unpublished) result on nodal sets when the geodesic flow of (M, g) 
is ergodic. Since there exist many expositions of quantum ergodic eigenfunctions, we only 
briefly recall the main facts and definitions and refer to \Z5\ IZ6j for further background. 

Quantum ergodicity concerns the semi-classical (large A) asymptotics of eigenfunctions 
in the case where the geodesic flow G* of (M, g) is ergodic. We recall that the geodesic 
flow is the Hamiltonian flow of the Hamiltonian H{x,C,) = |^|^ (the length squared) and 
that ergodicity means that the only G*-invariant subsets of the unit cosphere bundle S*M 
have either full Liouville measure or zero Liouville measure (Liouville measure is the natural 
measure on the level set = 1 induced by the symplectic volume measure of T*M). 

We will say that a sequence {fj^} of L^-normalized eigenfunctions is quantum ergodic if 

(34) (^^.., ^..) ^ j^^^^ ^Adfi, VA e ^°(M). 

Here, ^'^(M) denotes the space of pseudodifferential operators of order s, and dfi denotes 
Liouville measure on the unit cosphere bundle S*M of (M, g). More generally, we denote by 
dfir the (surface) Liouville measure on dB*M, defined by 

(35) dfir = ^ on dBlM. 

We also denote by a the canonical action 1-form of T*M. 

The main result is that there exists a subsequence {^j^} of eigenfunctions whose indices 
jk have counting density one for which Pj^{A) := {A(pj^,(pji^) — )• u{A) (where as above 
'^(^) ~ fj.{J*M) ^s•'M^^^^^ normalized Liouville average of a a)- Such a sequence of 

eigenfunctions is called a sequence of 'ergodic eigenfunctions'. The key quantities to study 
are the quantum variances 

(36) Va{X) := ^ E KM,^.) -^(^)r- 

The following result is the culmination of the results in |Sh.H [ZT| \CY\ IZZwj IGLj . 
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Theorem 4.4. Let {M,g) be a compact Riemannian manifold (possibly with boundary), and 
let {Aj,</9j} be the spectral data of its Laplacian A. Then the geodesic flow G* is ergodic on 
{S*M, (iyu) if and only if, for every A G \1/°(M), we have: 

(1) limA_oo^A(A) = 0. 

(2) (Ve)(35)limsup;,^^ ^^,>.fc:A,,A,,<A \{A^j,^k)\'^ < e 

Since all the terms in (1) are positive, no cancellation is possible, hence (1) is equivalent 
to the existence of a subset 5 C N of density one such that Qs ■= {d^k ■ k E S} has only cu 
as a weak* limit point. 

We now consider nodal sets of quantum ergodic eigenfunctions. The following result says 
that if we equip nodal sets with the measure Y2\V<fx.\dS , then nodal sets cannot 'scar', i.e. 

concentrate singularly as Aj — )■ oo. 

Proposition 1. Suppose that {v^a^} is a quantum ergodic sequence. Then any weak limit 
of {Yi\V(px \dS} must be absolutely continuous with respect to dV . 

We recall that, for any / G C2(M), 

(37) / ((A + A^)/) \^x\dV = [ flVifxldS. 

The identity for general / G C^(M) can be used to investigate the equidistribution of 
nodal sets equipped with the surface measure [Vv^'aM'S'. We denote the normalized measure 

Lemma 2. The weak * limits of the sequence {X~'^\^ ^\j\dS\z^^} of bounded positive measures 
are the same as the weak * limits of {\ipx.\\ (against f G C{M).) 



We let / G C (M) and multiply the identity (37) by A . We then integrate by parts to 
put A on /. This shows that for / G C^(M), we have 



f f\^x\dV = \-^ I f\Vx\dS + 0{X 

JM Jz^^ 



Letting / = 1, we see that the family of measures {A"^|Vv2Aj ^^(v^Aj )} is bounded. By 
uniform approximation of / G C{M) by elements of C^(M), we see that the weak* limit 
formula extends to C{M). 

Lemma 3. Suppose that {v^Aj} o quantum ergodic sequence. Then any weak limit of 
{\ip\^\dS} must be absolutely continuous with respect to dV . 

We recall that a sequence of measures /i„ converges weak * to if /^^ fd^in — ^ / fd^i for 
all continuous /. A basic fact about weak * convergence of measures is that / fd^n J fdfj, 
for all / G C(M) implies that fin{E) fi{E) for all sets E with fJ,{dE) = (Portmanteau 
theorem) . 

We also recall that a sequence of eigenfunctions is called quantum ergodic (in the base) if 
(38) I f^^^'^'"'' ^ V^Jj'"'- 
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In other words, (f\^lm the weak * topology, i.e. the vague topology on measures. We 
now prove Lemma [3j 

Proof. Suppose that |(y9Ajj.M^ dfi and assume that dfi = cdV + du where du is singular 
with respect to dV. Let S = supp u, and let a = = Let Te be the e-tube around 

S. Then 

lim [ \ipx \dV = cVol{T,) + = a + 0{e). 

fc— 5-00 Jj-^ 

But for any set C M, J^\vx^\dV < ^Vol{Vl)^ W\,?dV . Hence if Vol{dn) = 0, 
limsupj.,^ \ipx,\dV < Vol{n). Letting n = we get a + 0(e) < Vol{T,{J:)) = 0(e) 

since limfc_>oo J-j-^ Iv^Aj^. P*^^ = Vol (71) = 0(e). Letting e — )■ gives a contradiction. 

□ 

Of course, it is possible that the only weak* limit is zero. 

4.7. Weak* limits for quantum ergodic sequences. To our knowledge, the ques- 
tion whether the limit (4.4) holds / G L°° when (M, g) has ergodic geodesic flow has not been 
studied. It is equivalent to strengthening the Portmanteau statement to all measurable sets 
E, and is equivalent to the statement that {fj,.} — ^ 1 weakly in L^. We call such sequences 
L°° quantum ergodic on the base. The term 'on the base' refers to the fact that we only 
demand quantum ergodicity for the projections of the 'microlocal lifts' to the base M. For 
instance, the exponential eigenf unctions of flat tori are L°° quantum ergodic in this sense. 

Lemma 4.5. Suppose that {fj} is an quantum ergodic sequence. Then there exists 
e > so that \\<fj\\L'^>e>0 for all j . 

Proof. We argue by contradiction. If the conclusion were false, there would exist a subse- 
quence ipj^ — )■ strongly in L^, but with y^'j^.dV — )■ dV weakly in L^. The first assumption im- 
phes the existence of a subsequence (which we continue to denote by ipj^,) satisfying ^9^^. — )■ 
a.e. dV. But has the weak Banach-Saks property: any weakly convergent sequence in 
has a subsequence whose arithmetic means converge strongly (Szlenk's weak Banach-Saks 
theorem for L^). We choose such a subsequence for ipj^ and continue to denote it as ipj^. 
This subsequence has the properties that 

(1) fjk a.e. 

(2) ipN ■■= jf Y.k<N^h 1 strongly in L^. 

But TpNix) — )■ on the same set where '{>j^{x) — )■ 0, hence by (1) ipN a.s. This 
contradicts (2) and completes the proof. 

□ 

Combining with the above, we have 

Corollary 4.6. Suppose that {v^Aj } is o,n quantum ergodic sequence on the base. Then 
the conjectured Yau lower bound holds: 'H'^~'^{Z^^) > CgX for some Cg > 0. 

We also see that the limits in Proposition [T] are non-zero: 

Corollary 4.7. Suppose that {fx^} is an L°° quantum ergodic sequence on the base. Then 
there exists C > so that any weak limit of the sequence -^{Vipx^ldSlz^^ has mass > C > 0. 
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Of course, such an abstract functional analysis argument only serves a purpose if we can 
prove that eigenfunctions of A are quantum ergodic on the base in interesting cases. It 
is natural to conjecture that this condition holds on negatively curved manifolds, since the 
expected norm of a random wave is bounded below by a positive constant. The main 
problem is that L°°{M) is a non-separable Banach space. The standard quantum ergodicity 
arguments show that (when quantum ergodicity is valid), for any Borel set E there exists a 
subsequence iS^; of density one so that 



(39) 1™ . / = ^«^(^)- 

However, the non-separability of L°°(M) means that one cannot use the diagonalization 
argument of \L1\ ICVj to show that there exists a density one subsequence independent of 



E so that (39) holds. If quantum ergodicity fails, then zero-density subsequences of 
eigenfunctions would 'scar' along Cantor sets C of positive measure. That is, the mass 
J^ip'j^dV may tend to a larger value than Vol{C). 

Equidistributed sums of Gaussian beams and quantum ergodicity 

We briefly consider the question whether it is possible to have a quantum ergodic sequence 
of eigenfunctions for which Hv^jIIli — )• 0. 

First, we observe that there do exist sequences of quantum ergodic functions (not eigen- 
functions) with this property: Y.fJ? ^]w^X[a;^(„),x,(n)+i] ^ in L^{[0,l],dx) as long as 

M(n) = o{n). But its square is the probability measure j^^J2j=i ''''^X[xj(n),xj{n)+^ if 
the {xj{n)} are uniformly distribution in [0, 1] (w.r.t. dx), this tends weakly to dx. 

It is tempting to construct sequences of eigenfunctions with the same property: a Gaussian 
beam on the standard S'^ associated to a closed geodesic 7 (i.e. a rotate of Y^) is of 

height in a tube of radius y/X around 7. If we let M(N) = o{N^) and choose M(N) 
closed geodesies which are ^^^^^^ -separated, and become equidistributed in the space of 

closed geodesies, then (fN = / ^ ^^=f^ is an eigenfunction whose L-^-norm tends to 

zero like ^^M{N)N-i but whose norm is asymptotic to 1 and whose modulus square 
tends weak* to 1. More precisely, I^^P ~^ 1 weakly. To prove that Iv^Afp 1 

requires proving that j^f^ J^j^k'^-y^'^-y^ ~^ 0- '^^^ ~ M{N)^ terms which are 

exponentially outside the tube intersections T 1 (7j)nT 1 (7;^.). In the sum we may fix j = jo 
and multiply by M{N) . So we need then to show that X^fc^j^ | (K^^ , Y^ ) | — )■ 0. The geodesies 
are well-separated if the distance in the space of geodesies between them is > which 

means that the angle between 7^ and 7a; is at least this amount. When the angle is > e 
then the inner product | (Ki^, Y^) \ < -N~^ since the area of T 1 (7,) fl T 1 (7^) is bounded 
by this amount. For any e the sum over geodesies separated by e is 0{^M{N)N~^). The 
remaining number of terms is 0{e^M{N)). So if e = o(a/M(A^)) both terms tend to zero. 



4.8. Intersections of nodal sets of orthogonal eigenfunctions. A related question is 
whether nodal sets of orthogonal eigenfunctions of the same eigenvalue must intersect. Of 
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course, this question only arises when the eigenvalue has multiplicity > 1. A result of this 
kind was obtained by V. Gichev under a topological condition on M. 

Theorem 4.8. [Gi] Suppose that H^{M) = and that v5a,15V^a,2 o'^e orthogonal eigenfunc- 
tions with the same eigenvalue \^ . Then Z^^ ^ fl Z^^ 2^^- 

We briefly sketch the proof: Let A\ resp. A2 be the family of nodal domains of (/9a,i 
resp. ^\2- Each union Uive^j ^ covers M up to the nodal set of If the nodal sets 
do not intersect then the nodal set of v9a,2 is contained in Uvye^i instance; similarly 

if the indices are reversed. Hence the nodal sets have empty intersection if and only if 
Uw'eyli ^ UH/ey^l2 ^ covers M. Under this condition, Gichev constructs a closed 1-form 
which is not exact by showing that the incidence graph of the cover obtained from the union 
of the nodal domains of 93a,i and 99a,2 contains a cycle. He then considers a nodal domain V 
of v9a,i and a nodal domain V of (f\^2 which intersect. Let Q = dU fl V. Since Q fl dV 7^ 
there exists a smooth function f on M such that / = 1 in a neighborhood of Q and / = 
near dU\Q. Let 77 be the one form which equals df on U and on the complement of U. 
Clearly 77 is closed and it is verified in [Gi] that rj is not exact. 

Givech also proves that for S"^, if is a regular value of v9a,i then i^Z^,^ ^ fl Z^^ ^ ^ 2 for 
every orthogonal eigenfunction (y9A,2 with the same eigenvalue. The proof is simply to use 
Green's formula for a nodal domain for (px,i and note that the integral of v^a,2J^V^a,i equals 
zero on its boundary. 

A related observation is the curious identity of [SoZ], which holds for any {M,g): for any 
pair of eigenfunctions, 

(A^^-A^) [ ^aJ^a,MI^ = 2 / ifx,\Vipx,\dS. 

Hence for a pair of orthogonal eigenfunctions of the same eigenvalue, 

/ v^aJVv^aJc^'S' = 0. 

5. Norms and nodal sets 

Studies of nodal sets often involve dual studies of norms of eigenfunctions. In this 
section, we review a number of relatively recent results on norms, both in the global 
manifold M and for restrictions of eigenfunctions to submanifolds. 

5.1. Polterovich-Sodin on norms and nodal sets. Let ^(v^a) denote the collection of 
nodal domains of <^a- For A G A{ipx) let niA = max^ I^^aI- In |PS] the following is proved 
(see Corollary L7): 

Theorem 5.1. |PSj Let {M,g) be a C°° Riemannian surface. For every ip\ with ||v?a|| = 1; 

AaA 

Hence, for each a > 0, the number of nodal domains A of ipx where the maximal bound 
tha > aA^/^ is achieved in order of magnitude does not exceed kga~^. In particular, for fixed 
a, it remains bounded as A 00. 
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The proof uses a certain Bananch indicatrix, the Sogge bounds, and estimates on the 
inradius of nodal domains. For a continuous function u G C(]R), the generalized Banach 
indicatrix is defined by 

/ + 00 
n(c)/3(c,/)dc, 
-oo 

where for a regular value c G M of /, (3{c,f) is the number of connected components of 
/~^(c). In |PSj . the integral B{u,f)) is bounded from above through the L^-norms of the 
function / and A/. I.e.. in Theorem 1.3. For any / G J-a and any continuous function u on 

M, 

Biu,f)<k,\\uomf\\ + \\Af\\). 

The proof is roughly as follows: Let Pi be a point of Ai where th maximum is achieved. By 
the inradius bound jMan3j . there exists yU > so that the disc D{pj, j) (Z A^. One can then 
express (fx in D{pj, j) by the sum of a Green's integral and Poisson integral with respect 
to the Euclidean Dirichlet Green's function of a slightly smaller disc. In particular one may 
express fx{Pj) by such an integral. Apply Holder's inequality one gets 

JD{p,,r) 

Since the discs are disjoint one can sum in j and apply the Sogge bound to include the 
proof. Thus, the only fact one used about nodal domains was lower bound on the inradius. 

This result bears a curious comparison to the results of |STZj giving new constraints on 
(M, g) which are of maximal eigenfunction growth, i.e. possess eigenfunctions such that 
niA > CX2 for some sequence of eigenfunctions ipx^ with Xj — )■ 00. The result (building on 
older results of Sogge and the author) states that such a sequence can exist only if {M,g) 
possesses a 'pole' p for which the set of geodesic loops Cp based at p has positive measure in 
S*M (with respect to the natural spherical volume measure) and such that the first return 
map has a recurrence property. In fact, the only known surfaces where the bounds are 
achieved are surfaces of revolution, and in this case the first return map is the identity. It is 
quite plausible that if (M, g) has maximal eigenfunction growth, then the first return map 
must be the identity map on a set of positive measure in Cp. 

Combined with the Polterovich-Sodin result above, we see that such 'poles' p, when they 
exist, can only occur in a uniformly bounded number of nodal domains of a surface. It 
would be interesting to know if there can exist only a finite number of such points at all if 
one additionally assumes that the set of smoothly closed geodesies has measure zero. For 
instance,, in that case, there might be a unique pole in each of the finite number of possible 
nodal domains. This finitude problem would be useful in strengthening the condition on 
(M, g) of maximal eigenfunction growth. 

5.2. Norms of restrictions. A problem of current interest is to consider norms of 
restrictions of eigenfunctions to hypersurfaces or higher codimension submanifolds. For 
expository purposes we only consider geodesies on surfaces here. Following earlier work of 
A. Reznikov, Burq, Gerard and Tzvetkov |BGT] proved 

Theorem 5.2. |BGT] Suppose that {M,g) is a compact surface, then there exists Ao(e),C > 
so that, for any geodesic segment 7 of length and any eigenfunction ipx with X> Xq we 
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have 

(40) [ \^x\'ds<CX'^\\ipx\\^ 

Their estimate is sharp for the round sphere S*^ because of the highest weight spherical 
harmonics They also showed that for all geodesic segments 7 of unit length, 

The estimate is only known to be achieved when the geodesic is elliptic, and quite likely 
it can be improved if the geodesic is hyperbolic. A result in this direction is: 

Theorem 5.3. |SoZ2] Suppose that {M,g) is a compact surface of non-positive curvature. 
Then for all e, there exists Ao(e), C > so that, for any geodesic segment 7 of length and 
any eigenfunction ip\ with A > Ao(e), we have 

(41) / M'ds<Ce\'2\\^,\\^ 

A related result on norms is. 
Theorem 5.4. |SoZ3j Let {M,g) be a surface and assume that the set 

(42) r = {{x,0^ S*M : g\x, = {x, 0, some t > 0} 

of periodic points has Liouville measure zero in S*M. Then there is a subsequence of eigen- 
values Xj^ of density one so that 

(43) \\exjLHM) = o{Xy;). 

The results are based in part on a relatively new Kakeya-Nikodym maximal function 
estimate of Bourgain |Bourg| , as improved by Sogge |Sog2| . We believe that it can be 



improved the following phase space Kakeya-Nikodym theorem. Let Ts{'y) be the tube of 
radius 6 around a geodesic arc in M, and let x<5,7 be a smooth cutoff to a phase space tube 
of its lift to S*M. Then for all e, there exists 6{e) such that 



limsup— ^ sup / Iv^Apc^s < e. 
I\[X) f-< Ten Jt,,,w., 



We expect the sup occurs when 7 is the orbit of {x,C,). But then it is easy to estimate the 



right side and one should be able to get a quantitative improvement of Theorem 5.4 



5.3. Quantum ergodic restriction (QER) theorems. In this section we briefly review a 
recent series of results |TZ2t ITZ3i IDZl ICTZj on quantum ergodic restriction theorems. They 



are used in section ^ 10 to determine the limit distribution of intersections of nodal lines and 
geodesies on real analytic surfaces (in the complex domain). 

Let if C M be a hypersurface and consider the Cauchy data {^jIh, Xj^d^Lpjln) of eigen- 
functions along H; here di, is the normal derivative. We refer to ^PjIh as the Dirichlet data 
and to Xj^du^jln as the Neumann data. A QER (quantum ergodic restriction) theorem seeks 
to find limits of matrix elements of this data along H with respect to pseudo-differential op- 
erators Opnia) on H. The main idea is that S}jM, the set of unit covectors with footpoints 
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on H, is a cross-section to the geodesic flow and the first return map of the geodesic flow for 
S^M is ergodic. The Cauchy data should be the quantum analogue of such a cross section 
and therefore should be quantum ergodic on H. 

For applications to nodal sets and other problems, it is important to know if the Dirichlet 
data alone satisfies a QER theorem. The answer is obviously 'no' in general. For instance 
if (M, g) has an isometric involution and with a hypersurface H of fixed points, then any 
eigenfunction which is odd with respect to the involution vanishes on H. But in |TZ2l ITZ3j 
a sufficient condition is given for quantum ergodic restriction, which rules out this and 
more general situations. The symmetry condition is that geodesies emanating from the 'left 
side' of H have a different return map from geodesies on the 'right side' when the initial 
conditions are reflections of each other through TH. To take the simplest example of the 
circle, the restriction of sin kx to a point is never quantum ergodic but the full Cauchy data 
(cos fcx, sin fcx) of course satisfies cos^ kx + sin^ fcx = 1. In |CTZj it is proved that Cauchy 
data always satisfies QER for any hypersurface. This has implications for (at least complex) 
zeros of even or odd eigenfunctions along an axis of symmetry, e.g. for the case of Maass 
forms for the modular domain 5'L(2,Z)/'H^ (see ^10). 

To state the QER theorem, we introduce some notation. We put 

(44) T*^M = {{q,OeT;M, q E H}, T* H = {iq,r]) e T* H, q E H}. 
We further denote by tth : T^M — > T*H the restriction map, 

(45) TrH{x,0=^\TH. 

For any orientable (embedded) hypersurface H C M, there exists two unit normal co- 
vector fields i>± to H which span half ray bundles N± = M+iz-t C N*H. Infinitesimally, they 
define two 'sides' of H, indeed they are the two components of T^M\T*H. We use Fermi 
normal coordinates (s, ?/„) along H with s E H and with x = exp^ ?/„z/ and let a, rjn denote 
the dual symplectic coordinates. For (s,cr) E B*H (the co-ball bundle), there exist two unit 
covectors C,±{s,a) E S*M such that |^±(s,cr)| = 1 and ^\tsH = ct- In the above orthogonal 
decomposition, they are given by 



(46) ^±{s,a) = a±y/l-\a\^u4s). 
We define the refiection involution through T*H by 

(47) th-.T^hM ^T*hM, rH{s,^^i±{s,a)) = {s,^^i^{s,a)), G M+. 

Its fixed point set is T*H. 

We denote by G* the homogeneous geodesic fiow of {M,g), i.e. Hamiltonian flow on 
T*M — generated by \^\g. We define the first return time T{s,^) on S*^M by, 

(48) T(s,O = inf{t>0:G*(s,OG^HM, {s,i) E S*hM)}. 

By definition T{s,^) = +oo if the trajectory through (s,^) fails to return to H. Inductively, 
we define the jth return time T^^\s,$,) to S*^M and the jth return map when the return 
times are finite. 

We define the first return map on the same domain by 

(49) $:^^M^^^M, $(s,e) = G^(^'«Hs,0 

When G* is ergodic, $ is defined almost everywhere and is also ergodic with respect to 
Liouville measure ^il^h on S*^M. 
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Definition: We say that H has a positive measure of microlocal reflection symmetry if 

Vjyo / 
Otherwise we say that H is asymmetric with respect to the geodesic flow. 

The QER theorem we state below holds for both poly-homogeneous (Kohn-Nirenberg) 
pseudo-differential operators as in |HoI-IVj and also for semi-classical pseudo-differential 
operators on H |Zwj with essentially the same proof. To avoid confusion between pseu- 
dodifferential operators on the ambient manifold M and those on H, we denote the latter 
by Oph{o.) where a G S^i{T*H). By Kohn-Nirenberg pseudo-differential operators we mean 
operators with classical poly-homogeneous symbols a{s,a) G C°°(T*H), 

oo 

a(s, a) ~ a_fc(s, a), (a_fc positive homogeneous of order — k) 

k=0 

as \a\ — )■ oo on T*if as in |HoI-IVj . By semi-classical pseudo-differential operators we mean 
/i-quantizations of semi-classical symbols a G S^'^{T*H x [0,ho]) of the form 

oo 

ah{s,a) ~ 

k=0 

as in [Z^lHZllTZ] . 

We further introduce the zeroth order homogeneous function 

II I |2 

(50) j{s,yn,a,vn)= , L^"' , , = (1 - ^)^ {r' = \a\' + \rjn\') 
on T^M and also denote by 

(51) 7^.^ = (i-|a|2)i 

its restriction to S*^M = {r = 1}. 

For homogeneous pseudo-differential operators, the QER theorem is as follows: 

Theorem 5.5. |TZ|, ITZ2| IDZ] Let {M,g) be a compact manifold with ergodic geodesic flow, 
and let H d M he a hypersurface. Let ^p\j;j = 1,2, ... denote the -normalized eigenfunc- 
tions of Ag. If H has a zero measure of microlocal symmetry, then there exists a density-one 
subset SofN such that for Aq > and a(s, a) G S^i(T*H) 



where 



lim {OpH{a)-fH^x,,lH^x,)L^H) = (^{a) 



V0L{O IVl) Jb*H 



Alternatively, one can write uj{a) = ^^k^^'M) Is* m ^'^^^'''^^^0)d'fJ'L,H{0- Note that ao{s,a) 
is bounded but is not defined for a = 0, hence ao(s, 7r//((^)) is not defined for ^ G N*H 
if ao(s,(j) is homogeneous of order zero on T*H. The analogous result for semi-classical 
pseudo-differential operators is: 
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Theorem 5.6. \TZ\ ITZ2| IDZ] Let {M,g) be a compact manifold with ergodic geodesic flow, 
and let H G M be a hypersurface. If H has a zero measure of microlocal symmetry, then 
there exists a density-one subset SofN such that for a G S'°''^(T*i!f x [0, ho)), 

lim {OphAa)^HVhplHVh,) L^H) = w(a), 

where 

/ ao(s,cr)7B.^//(s,(x) dso?(J. 

' B*H 

Examples of asymmetric curves on surfaces in the case where (M, g) is a finite area hy- 
perbohc surface are the following: 

• if is a geodesic circle; 

• if is a closed horocycle of radius r < inj{M,g), the injectivity radius. 

• if is a generic closed geodesic or an arc of a generic non-closed geodesic. 

6. Critical points 



vol{S*M) 



In this section, we briefly discuss some analogues of (16) and (24) for critical points on 
surfaces. To be sure, it is not hard to generate many identities; the main problem is to derive 
information from them. 

We denote the gradient of a function ip by Vip and its Hessian by V^v^ := Vdip, where V 
is the Riemannian connection. We also denote the area form by dA and the scalar curvature 
by K. The results are based on unpublished work in progress of the author. It is often said 
that measuring critical point sets and values is much more difficult than measuring nodal 
sets, and in a sense the identities reflect this difficulty, and we immediately see one difficulty 
in that the identities become signed: 

Proposition 6.1. Suppose that {M,g) is a Riemannian surface, and that if is a Morse 
eigenfunction with {A + X'^)(p = 0. Let V e C'^{M). Then 

2vr Ep:.,(p)=o ^^Wdet VV(p)) V{p) = 2X' \^\^dA + 2 KVdA 

(52) 

-J,,iAV)log\V^\'dA. 

Here, sign(det V'^ip{p)) = 1 if p is a local maximum or minimum and = —1 if p is a saddle 
point. When V = 1, the identity reduces to the Gauss-Bonnet theorem J KdA = 2ttx{M) 
and the Hopf index formula x(^) = ^x\/ip{x)=o^^S^i^^^^'^'^iP))- -^^ ^^^^ indicates, the 
main problem with applying the identity to counting critical points is that the left side is an 



alternating sum over critical points rather than a positive sum. In Dong a related identity 
using I V<y5p + A^v9^ produced a sum of constant sign over the singular points of ip, but singular 
points are always saddle points of index —1 and hence of constant sign. Note that under 
the Morse assumption, log iVy^l, iVy?!""^ G L^{M, dA), so that the right side is a well defined 
measure integrated against V . 

We now make some interesting choices of V . As mentioned above, (weighted) counting 
of critical values should be simpler than weighted counting of critical points. Hence we put 
V = f{ip) for smooth /. This choice does give cancellation of the 'bad factor' |Vv5|~^ and 
(using that A/(^) = r(^)|V<^p - /'(<^)AV) we get 
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Corollary 4. With the assumptions of Proposition 6.1, if f e C^(]R), then 
(53) 

- /m(/"(^)I V^P - f(^)AV)) log \V^\'dA. 

Of course, this still has the defect that the left side is an oscillating sum, and the factor 
f{^p) in the sum damps out the critical points in regions of exponential decay. To illustrate, 
if f{x) = X we get 

(54) 27rEp;,^(p)=oSign(detVV(p))^(p) = 2 K^dA + f^^ ^\og\V^\^dA. 

To highlight the sign issue, we break up the sum into the sub-sum over maxima/minima 
and the sub-sum over saddle points, denoting the set of local maxima (resp. minima) by max 
(resp. min) and the set of saddle points by Sad. Of course we have #(max U min) — ^Sad = 



x(M). Then (53) is equivalent to 



(55) 

2^Ep„umin)/(v^(p)) = 2vrEpe5a.)/(v^(p)) + 2X' ^ f i^)dA + 2 J^^ K f i^)dA 

- /,,(/"(¥.) I Vy^p - r(¥.)AV)) log I V^prfA. 

We write logr = log_,_r — log_ r where log_,_r = max{logr, 0}. We note that on any 
compact Riemannian manifold, log^ |Vv?p = O(logA) uniformly in x as A — )■ oo while 
log_ |V(/5p can be quite complicated to estimate. When f = we get, 

(56) 

+4A2 + 2 X,,(AV' - I V</^n log \Vif\^dA + 2 K<f^dA. 
Assuming cp is a. Morse eigenfunction, this implies 

(57) Zpen.„^'iP) < EpeSaMpf + OiXHogX). 

To get rid of the signs in the sum, we could choose V = W det V^v', where the determinant 
is defined by the metric. Since (sign det V^v^) det V^v' = | det V^v^| we obtain 

(58) 

2^Ep:M.)=o|detVV(p)|W^(p)) = 2A^ ^ v(H^det^v^,).v. ^^ + 2 KW det V^^dA 

- f^i^W det VV) log \V^\^dA. 
But the first term appears to be difficult to estimate. 

7. Analytic continuation of eigenfunctions for real analytic {M,g) 



We now take up the theme mentioned in the introduction of analytically continuing eigen- 
functions on real analytic (M, g) to the complex domain. In the next sections we apply the 
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analytic continuation to the study of nodal of eigenfunctions in the real analytic case. For 
background we refer to [LSll [GSH [GS2l iGLSl [Z8] . 

A real analytic manifold M always possesses a unique complexification M<c generalizing the 
complexification of as C™. The complexification is an open complex manifold in which 
M embeds l : M ^ Mq as a totally real submanifold (Bruhat- Whitney). As examples, we 
have: 

• M = W^/Z"" is Mc = C'^/Z"'. 

• The unit sphere S"' defined by + ■ ■ ■ + x'^^i = 1 in ]R"+^ is complexified as the 
complex quadric = {{zi, . . . , z^) G C""^"*^ : z'l + ■ ■ ■ + z^^-^ = 1}. 

• The hyperboloid model of hyperbolic space is the hypersurface in ]R"+^ defined by 

W = {xl + ---xl- X^+i = -1, Xn > 0}. 

Then, 

= {{z„ G C"+i :zl + ...zl- zl^, = -1}. 

• Any real algebraic subvariety of M"* has a similar complexification. 

• Any Lie group G (or symmetric space) admits a complexification Gc. 

The Riemannian metric determines a special kind of distance function on known as a 
Grauert tube function. It is the plurisubharmonic function ^ = on defined as the 
unique solution of the Monge- Ampere equation 

(ddy^)"' = SM^,dVg, t*{iddp) = g. 

Here, SM^,dVg is the delta-function on the real M with respect to the volume form dVg, 
i.e. / — >■ Jj^j fdVg. In fact, it is observed in |GSH IGLS] that the Grauert tube function 

is obtained from the distance function by setting y/p{C) = ^\/r^(C, C) where r'^{x,y) is the 
squared distance function in a neighborhood of the diagonal in M x M. 

One defines the Grauert tubes M^- = {( E Mq : y^(C) < t}- There exists a maximal tq 
for which is well defined, known as the Grauert tube radius. For r < tq, Mj. is a strictly 
pseudo-convex domain in Mq. 

The complexified exponential map (x, ^) — )■ exp^i^ defines a diffeomorphism from B*M to 
Mr and pulls back ^/p to \^\g. The one-complex dimensional null foliation of dd^, known 
as the 'Monge- Ampere' or Riemann foliation, are the complex curves t + — )■ T'y{t), where 
7 is a geodesic, where r > and where T'jit) denotes multiplication of the tangent vector to 
7 by r. We refer to |LSlt IGLSl [Z8] for further discussion. 

7.1. Poisson operator and analytic Continuation of eigenfunctions. The half- wave 
group of (M, g) is the unitary group U (t) = e'*^ generated by the square root of the positive 
Laplacian. Its Schwartz kernel is a distribution on M x M x M with the eigenfunction 
expansion 

oo 

(59) U{t,x,y) = Y,e''''fjix)fj{y)- 

j=0 

By the Poisson operator we mean the analytic continuation off/ (t) to positive imaginary 
time, 

(60) e""^ = U{iT). 
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The eigenfunction expansion then converges absolutely to a real analytic function on IR+ x 
M X M. 

Let A{t) denote the operator of analytic continuation of a function on M to the Grauert 
tube Mr- Since 

(61) Uc{tT)^x = c-^Va, 
it is simple to see that 

(62) A{t) = f/c(2r)e"^ 

where Uc{iT, (, y) is the analytic continuation of the Poisson kernel in x to M^-. In terms of 
the eigenfunction expansion, one has 

oo 

(63) U{iT,(,y) = J2^'^''^J^0vM, (Cz/)eM,xM. 

j=0 

This is a very useful observation because Uc{iT)e'^^ is a Fourier integral operator with 
complex phase and can be related to the geodesic flow. The analytic continuability of the 
Poisson operator to implies that every eigenfunction analytically continues to the same 
Grauert tube. 

7.2. Analytic continuation of the Poisson wave group. The analytic continuation of 
the Possion-wave kernel to M^- in the x variable is discussed in detail in |Z8j and ultimately 
derives from the analysis by Hadamard of his parametrix construction. We only briefly 
discuss it here and refer to |Z8j for further details. In the case of Euclidean M" and its wave 
kernel U{t,x,y) = e**l^le*^^'^~^^(i(^ which analytically continues to t + irX = x + ip G 
C+ X C" as the integral 

Ucit + iT,x + ip, y)= [ e^(*+^")l«le*<«'^+^P-^>rfe 



The integral clearly converges absolutely for \p\ < r. 

Exact formulae of this kind exist for 5™ and H™. For a general real analytic Riemannian 
manifold, there exists an oscillatry integral expression for the wave kernel of the form, 

(64) Uit,x,y)= [ e**l«l««e'<«'""P«'(")U(t,a;,?/,Ot^e 

Jt*m 

where A{t,x,y,^) is a poly homogeneous amplitude of order 0. The holomorphic extension 
of rt64| to the Grauert tube \C\ < r in a; at time t = ir then has the form 



(65) Uc{tT,C,y)= e-^l«l«^e^<«'^-P^"'('^)U(t,C,l/,Orfe {( = x + tp). 



T* 
y 



7.3. Analytic continuation of eigenfunctions. Thus, a function / G C°°{M) has a 
holomorphic extension to the closed tube ^/p{C) < r if and only if / G Dom{e^^), where 
gTv^ jg ^YiQ backwards 'heat operator' generated by (rather than A). That is, / = 
Yl^=o ^n^Xn admits an analytic continuation to the open Grauert tube if and only if / 
is in the domain of e"^^, i.e. if |a„pe^'^^" < oo. Indeed, the analytic continuation is 
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f/c(^r)e^^/. The subtlety is in the nature of the restriction to the boundary of the maximal 
Grauert tube. 

This result generalizes one of the classical Paley- Wiener theorems to real analytic Rie- 
mannian manifolds |BoutlGS2] . In the simplest case of M = S*^, / ~ ^^ez ctnC*"^ G C"^(S'^) is 
the restriction of a holomorphic function F ~ X^nez ^^z" on the annulus = {\ log |2;|| < r} 
and with F G L'^{dS\) if and only if \ f{f^)\'^ e^'"'^ < oo. The case of is more compli- 
cated since it is non-compact. We are mainly concerned with compact manifolds and so the 
complications are not very relevant here. But we recall that one of the classical Paley-Wiener 
theorems states that a real analytic function / on M" is the restriction of a holomorphic func- 
tion on the closed tube |ImC| < r which satisfies \F{x + i^)\'^dx < C for ^ < r if and 
only if /Vl^"^^l G L2(M"). 

Let us consider examples of holomorphic continuations of eigenf unctions: 

• On the flat torus W^/II^, the real eigenfunctions are cos(fc, x), sin(A;, x) with k G 
2tt7/^. The complexified torus is C^/Z™ and the complexified eigenfunctions are 
cos{k, , sm{k, () with ( = x + i^. 

• On the unit sphere S"^, eigenfunctions are restrictions of homogeneous harmonic 
functions on M™"*"^. The latter extend holomorphically to holomorphic harmonic 
polynomials on C™"*^^ and restrict to holomorphic function on S^. 

• On H™, one may use the hyperbolic plane waves e(*'^+^)^^'''\ where {z, b) is the (signed) 
hyperbolic distance of the horocycle passing through z and b to 0. They may be 
holomorphically extended to the maximal tube of radius 7r/4. 

• On compact hyperbolic quotients H™/r, eigenfunctions can be then represented by 
Helgason's generalized Poisson integral formula jH] . 

^^{z) = [ e'^'^+'^^'''^dTx{b). 
Jb 

Here, z E D (the unit disc), B = dD, and dTx G V'{B) is the boundary value of v^a, 
taken in a weak sense along circles centered at the origin 0. To analytically continue 
ipx it suffices to analytically continue {z,b). Writing the latter as (C,^), we have: 

(66) ^^(C) = / e(^^+i)<^-'')rfT,(6). 

JB 

7.4. Complexified spectral projections. The next step is to holomorphically extend the 
spectral projectors dU.]f)^x]{x-,y) = X]j ■^j)fji^)fjiy) of The complexified diagonal 
spectral projections measure is defined by 

(67) dxUf,^x]{C,C) = J2^{>^-^MiO\'- 

j 

Henceforth, we generally omit the superscript and write the kernel as H^ ;^j(^, Q. This kernel 

is not a tempered distribution due to the exponential growth of |v5j'(C)P- Since many as- 
ymptotic techniques assume spectral functions are of polynomial growth, we simultaneously 
consider the damped spectral projections measure 

(68) dxP[o,x]iCX) = E^(^ - ^3)e-'^"'\vfiO\\ 

3 
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which is a temperate distribution as long as ^/f>{C) ^ When we set r = ^/f>{C) we omit 
the r and put 

(69) c^,P[o,A](c,c) = E'^(^ - 1^5^(0 r- 

The integral of the spectral measure over an interval / gives 

j-.Xj&I 



Its complexification gives the kernel (118) along the diagonal, 

(70) MCC)= E 



and the integral of (68) gives its temperate version 

(71) ^'/(c,c)= E e"'^'^' 1^5^(0 1', 

or in the crucial case of r = y/p{(), 

(72) Piic,c)= E e-^^(^)'^ 1^,^(0 r, 

7.5. Poisson operator as a complex Fourier integral operator. The damped spectral 
projection measure d\ P[QX]i(X) (68) is dual under the real Fourier transform in the t 
variable to the restriction 

(73) Uit + 2tT,C,C) = 

j 

to the anti-diagonal of the mixed Poisson-wave group. The adjoint of the Poisson kernel 
U{iT,x,y) also admits an anti-holomorphic extension in the y variable. The sum (73) are 
the diagonal values of the complexified wave kernel 

U{t + 2tTX,C) = Jj,,U{t + iTX,y)E{iT,yX')dVg{x) 

(74) 

^^^.,(-2.+.t)A,^C(^)^C(^,). 

We obtain ( [74| by orthogonality of the real eigenfunctions on M. 

Since U{t + 2iT, (, y) takes its values in the CR holomorphic functions on dMr, we consider 
the Sobolev spaces O*"' ^(dM^-) of CR holomorphic functions on the boundaries of the 
strictly pseudo-convex domains M^, i.e. 

where Wg is the sth Sobolev space and where 0{dM^) is the space of boundary values of 
holomorphic functions. The inner product on 0^{dMr) is with respect to the Liouville 
measure 

(75) dfir = {idd^)""-^ A d^^. 
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We then regard U{t + iT,(,y) as the kernel of an operator from L'^{M) — )■ 0^{dMr). It 
equals its composition 11,- o U(t + it) with the Szegoprojector 

: L\dMr) -> 0\dMr) 

for the tube Mr, i.e. the orthogonal projection onto boundary values of holomorphic func- 
tions in the tube. 

This is a useful expression for the complexified wave kernel, because is a complex 
Fourier integral operator with a small wave front relation. More precisely, the real points 
of its canonical relation form the graph of the identity map on the symplectic one 
C T*dMj. spanned by the real one-form df^p, i.e. 

(76) = {(C;rf/V(C)), C e r > 0} C T\dMr). 

We note that for each r, there exists a symplectic equivalence ~ T*M by the map 
{(,rd'^p{Q) — )■ {E^^{Q,ra), where a = ^ ■ dx is the action form (cf. |GS2j ). 
The following result was first stated by Boutet de Monvel (for more details, see also 

Theorem 7.1. |Boul IGS2j n, o U{ie) : L'^{M) 0{dM^) is a complex Fourier integral 
operator of order — associated to the canonical relation 

r = {(y,r/,i,(2/,r7)}cT*MxS,. 

Moreover, for any s, 

n, o U{ie) : W'{M) O'+'^idM,) 
is a continuous isomorphism. 

In |Z8] we give the following sharpening of the sup norm estimates of |Boul IGLS] : 

Proposition 7.2. Suppose {M,g) is real analytic. Then 

sup \v^,{0\ < CX'^e^\ sup l^^l < CA^e-^ 

The proof follows easily from the fact that the complexified Poisson kernel is a complex 
Fourier integral operator of finite order. The estimates can be improved further. 

7.6. Maximal plurisubharmonic functions and growth of (p'^. In |Z8j . we discussed 
analogues in the setting of Gruaert tubes for the basic notions of pluripotential theory on 
domains in C". Of relevance here is that the Grauert tube function y/p is the analogue of the 
pluri-complex Green's function. We recall that the maximal PSH function (or pluri-complex 
Green's function) relative to a subset E G Q is defined by 

^^(0 = sup{n(z) : u G PSH{n),u\E < 0,u\9n < !}• 
On a real analytic Riemannian manifold, the natural analogue of is the space 

= {p= ^ aj(fx,, ai, . . . , aAr(A) G M} 

j:Xj<X 
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spanned by eigenfunctions with frequencies < A. Rather than using the sup norm, it is 
convenient to work with based norms than sup norms, and so we define 

7V(A) 

= {V= aj<^A,, ||p||L2(Af) = kiT = !}• 

r.\j<\ j=i 

We define the A-Siciak extremal function by 

^l,{z) = sup{|^(^)|i/^: ij e -Hx; Mm ^ 1}, 

and the extremal function by 

<!>m{z) = sup<l>^,^(z). 

A 

The extremal PSH function is defined by 

VgiCr) = snp{u{z) : u G PSH{Mr),u\M < 0,M|aAf. < r}. 
In |Z8] we proved that Vg = ^/p and that 
(77) $M = Vg. 



The proof is based on the properties of (70). By using a Bernstein- Walsh inequality 



N{x) - <^ucy 

it is not hard to show that 

(78) ^m{z) = lim |logn[o,A(C,C)- 

A— >oo A 

To evaluate the logarithm, one can show that the kernel is essentially e'^^ times the tem- 
perate projection defined by the Poisson operator, 

(79) P[o,A](c,c)= E e-2v^(^)^^i^5^(c)p. 

j:A,e[0,A] 



The equality (77) follows from the fact that limA-5>oo x ^oS-^[o,a](C) C) = 0- 



We now return to nodal sets, where we will see the same extremal functions arise. 
8. Counting nodal lines which touch the boundary in analytic plane 

DOMAINS 

It is often possible to obtain more refined results on nodal sets by studying their inter- 
sections with some fixed (and often special) hypersurface. This has been most successful in 
dimension two. In this section, we review the results of |TZj giving upper bounds on the 
number of intersections of the nodal set with the boundary of an analytic (or more generally 
piecewise analytic) plane domain. One may expect that the results of this section can also be 
generalized to higher dimensions by measuring codimension two nodal hypersurface volumes 
within the boundary. 

Thus we would like to count the number of nodal lines (i.e. components of the nodal set) 
which touch the boundary. Here we assume that is a regular value so that components of 
the nodal set are either loops in the interior (closed nodal loops) or curves which touch the 
boundary in two points (open nodal lines). It is known that for generic piecewise analytic 
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plane domains, zero is a regular value of all the eigenfunctions ip\.^ i.e. Vv^a^ 7^ on 
^vx- [U]? '^^ then call the nodal set regular. Since the boundary lies in the nodal set for 
Dirichlet boundary conditions, we remove it from the nodal set before counting components. 
Henceforth, the number of components of the nodal set in the Dirichlet case means the 
number of components of Z^p^ \^VL. 

In the following, and henceforth, > denotes a positive constant depending only on 
the domain VL. 

Theorem 8.1. Let Q be a piecewise analytic domain and let ngn(Aj) be the number of 
components of the nodal set of the jth Neumann or Dirichlet eigenfunction which intersect 
dQ. Then there exists such that ngQ{Xj) < C^Xj. 

By a piecewise analytic domain C M^, we mean a compact domain with piecewise 
analytic boundary, i.e. dQ is a union of a finite number of piecewise analytic curves which 
intersect only at their common endpoints. Such domains are often studied as archtypes 
of domains with ergodic billiards and quantum chaotic eigenfunctions, in particular the 
Bunimovich stadium or Sinai billiard. Their nodal sets have been the subject of a number 
of numerical studies (e.g. |BGSt IFGS] ). 

In general, there does not exist a non-trivial lower bound for the number of components 
touching the boundary. E.g. in a disc, the zero sets of the eigenfunctions are unions of circles 
concentric with the origin and spokes emanating from the center. Only the spokes intersect 
the boundary and their number reflects the angular momentum rather than the eigenvalue 
of the eigenfunction. But we conjecture that for piecewise analytic domains with ergodic 
billiards, the the number of complex zeros of (p'^ ldUc is bounded below by CnXj. We discuss 



work in progress on this conjecture in ^ 10 



In comparison to the order 0{Xj) of the number of boundary nodal points, the total number 
of connected components of Z^^ has the upper bound 0(A|) by the Courant nodal domain 
theorem. It is not known in general whether the Courant upper bound is achieved, but we 
expect that it is often achieved in order of magnitude. In [NSj it is proved that the average 
number of nodal components of a random spherical harmonic is of order of magnitude A^. 
Thus, the number of components touching the boundary is one order of magnitude below 
the total number of components. 

8.1. Boundary critical points. The article |TZj also contains a similar estimate on the 
number of critical points of ip\. which occur on the boundary. We denote the boundary 
critical set by 

C^,^={qEdn:id^,^)iq)=0}. 

In the case of Neumann eigenfunctions, q G C^^^ d^ipxjlanil)) = since the 

normal derivative automatically equals zero on the boundary, while in the Dirichlet case 
q G C^^, -^^^ dyip\.{q) = since the boundary is a level set. 

We observe that radial eigenfunctions on the disc are constant on the boundary; thus, 
boundary critical point sets need not be isolated. We therefore impose a non-degeneracy 
condition on the tangential derivative dt{fx^ \dn) to ensure that its zeros are isolated and can 
be counted. We say that the Neumann problem for a bounded domain has the asymptotic 
Schiffer property if there exists C > such that, for all Neumann eigenfunctions ip\. with 
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sufficiently large Xj, 



(80) " " > e 

Here, dt is the unit tangential derivative, and the norms refer to the restrictions of the 
eigenfunction to dQ. 

Theorem 8.2. Let f2 C 6e piecewise real analytic. Suppose that ip\ \dn satisfies the 
asymptotic Schiffer condition (80) in the Neumann case. Then the number of n^^j^^{\j) = 



4^Cip^ of critical points of a Neumann or Dirichlet eigenfunction ipx . which lie on dQ satisfies 



n 



cri 



ti^j) < CqXj for some Cq > 



In the case of Dirichlet eigenfunctions, endpoints of open nodal lines are always boundary 
critical points, since they must be singular points ofipx^. Hence, an upper bound for n^^^^{Xj) 
also gives an upper bound for the number of open nodal lines. 

Corollary 8.3. Suppose that Q (Z M.'^ is a piecewise real analytic plane domain. Let 
ng^{Xj) be the number of open nodal lines of the jth Dirichlet eigenfunction, i.e. connected 
components of {ip\^ = 0} C fl° whose closure intersects dfl. Then there exists Co > such 
that nQ^{Xj) < CnXj. 

There does not exist a non-trivial lower bound on the number of interior critical points 
[JNl. 



8.2. Proof by analytic continuation. For the Neumann problem, the boundary nodal 
points are the same as the zeros of the boundary values fx^lan of the eigenfunctions. The 
number of boundary nodal points is thus twice the number of open nodal lines. Hence in 



the Neumann case, Theorem 8.1 follows from: 



Theorem 8.4. Suppose that Q C M."^ is a piecewise real analytic plane domain. Then the 
number n{Xj) = i^Z^^ fl dQ of zeros of the boundary values ^px^lan of the jth Neumann 
eigenfunction satisfies n{Xj) < C^Xj, for some > 0. 



This is a more precise version of Theorem 8.1 since it does not assume that is a regular 



value. In keeping with the theme of this survey, we prove Theorem 8.4 by analytically 
continuing the boundary values of the eigenfunctions and counting complex zeros and critical 
points of analytic continuations of Cauchy data of eigenfunctions. When dQ G C^, the 
eigenfunctions can be holomorphically continued to an open tube domain in projecting 
over an open neighborhood W in of fl which is independent of the eigenvalue. We denote 
by flc C the points ( = x + i^ G with x G fi. Then (p\j{x) extends to a holomorphic 
function v^^^, (C) where x E W and where |^| < eo for some eo > 0. 

Assuming dfl real analytic, we define the (interior) complex nodal set by 

and the (interior) complex critical point set by 
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Theorem 8.5. Suppose that dM? is a piecewise real analytic plane domain, and denote 
by {dQ)c the union of the complexifications of its real analytic boundary components. 

(1) Let n(\j,d^lc) = ^Z^^"^ be the number of complex zeros on the complex boundary. 

Then there exists a constant Cn > independent of the radius of {dQ)c such that 



(2) Suppose that the Neumann eigenf unctions satisfy (80) and let n , dflc) = #C^^*^. 



Then there exists Cq > independent of the radiusoj {dVL)c such that n ^j.j^^{Xj , dilc) < 

CqXj. 

The theorems on real nodal lines and critical points follow from the fact that real zeros 
and critical points are also complex zeros and critical points, hence 

(81) n{Xj) < n{Xj,dQc)', '^crit(-^i) - '^crit(-^i' "^^c)- 

All of the results are sharp, and are already obtained for certain sequences of eigenfunctions 



on a disc (see ^4.3). If the condition (80) is not satisfied, the boundary value of (pxj must 
equal a constant Cj modulo an error of the form o{e~'"^^). We conjecture that this forces 
the boundary values to be constant. 



The method of proof of Theorem |8.5| generalizes from dQ to a rather large class of real 
analytic curves C G Q, even when dQ is not real analytic. Let us call a real analytic curve 
C a good curve if there exists a constant a > so that for all Xj sufficiently large, 

(82) WUw<e'^^.. 

II'^aJ|l2(c) 

Here, the norms refer to the restrictions of the eigenfunction to C and to dQ. The 
following result deals with the case where C C dQ is an interior real-analytic curve. The 
real curve C may then be holomorphically continued to a complex curve Cc C obtained 
by analytically continuing a real analytic parametrization of C. 

Theorem 8.6. Suppose that Q (Z M."^ is a C°° plane domain, and let C G Q be a good 



interior real analytic curve in the sense of (82). Let n{Xj,C) = i^Z^^ f] C be the number 



of intersection points of the nodal set of the j-th Neumann (or Dirichlet) eigenfunction with 
C. Then there exists Aq^q, > depending only on C,Q such that n{Xj,C) < Ac^nXj. 

A recent paper of J. Jung shows that many natural curves in the hyperbolic plane are 
'good' [JJ]. 



8.3. Application to Pleijel's conjecture. We also note an interesting application due to 



I. Polterovich |Poj of Theorem 8.1 to an old conjecture of A. Pleijel regarding Courant's 
nodal domain theorem, which says that the number of nodal domains (components of 
fl\Ztp^^) of the fcth eigenfunction satisfies Uk < k. Pleijel [P] improved this result for Dirichlet 
eigefunctions of plane domains: For any plane domain with Dirichlet boundary conditions, 
limsup^_5.o^ ^ < ^ — 0.691..., where ji is the first zero of the Jq Bessel function. He 
conjectured that the same result should be true for a free membrane, i.e. for Neumann 
boundary conditions. This was recently proved in the real analytic case by I. Polterovich 
[Poj . His argument is roughly the following: Pleijel's original argument applies to all nodal 
domains which do not touch the boundary, since the eigenfunction is a Dirichlet eigenfunction 
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in such a nodal domain. The argument does not apply to nodal domains which touch the 



boundary, but by Theorem 8.1 the number of such domains is negligible for the Pleijel bound. 



9. EQUIDISTRIBUTION of complex nodal sets of real ERGODIC EIGENFUNCTIONS 
ON ANALYTIC {M,g) WITH ERGODIC GEODESIC FLOW 

We now consider global results when hypotheses are made on the dynamics of the geodesic 
flow. Use of the global wave operator brings into play the relation between the geodesic 
flow and the complexified eigenfunctions, and this allows one to prove gobal results on nodal 
hypersurfaces that reflect the dynamics of the geodesic flow. In some cases, one can determine 
not just the volume, but the limit distribution of complex nodal hypersurfaces. Since we 
have discussed this result elsewhere |Z6j we only briefly review it here. 

The complex nodal hypersurface of an eigenfunction is defined by 

(83) Z^e = {C e B^M : y,^(C) = 0}. 

There exists a natural current of integration over the nodal hypersurface in any ball bundle 
B*M with e < eo , given by 

(84) ([Z^e],^) = ^ /" 9aiog|^^rA^= /" ^, ^GP(-i'™-i)(i?:M). 



B*M 



In the second equality we used the Poincare-Lelong formula. The notation p(™'^^'™^^)(i?*M) 
stands for smooth test (m — 1, m — l)-forms with support in B*M. 

The nodal hypersurface Z c also carries a natural volume form \Z c \ as a complex hyper- 



m — 1 



surface in a Kahler manifold. By Wirtinger's formula, it equals the restriction of to 
Z,„c. Hence, one can regard Z,„c as defining the measure 

(85) {\Z^c\,^) = j^ ^(^^^31)!' V^C{B:M). 

We prefer to state results in terms of the current [Z^c] since it carries more information. 

Theorem 9.1. Let {M,g) be real analytic, and let {^pj^} denote a quantum ergodic sequence 
of eigenfunctions of its Laplacian A. Let {B*^^M, J) be the maximal Grauert tube around M 
with complex structure Jg adapted to g. Let e < e^. Then: 

-^[Z c ] -ddy/p weakly m V'^^'^\B*M), 
in the sense that, for any continuous test form ip G T)^'^'^''^^^\B* M) , we have 

^3k JZ V ^ Jb*M 



3k 



Equivalently, for any ip G C{B*M), 



hk Jz c {m-iy. TT Jb*m (m - 1^ 
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Corollary 9.2. Let {M,g) be a real analytic with ergodic geodesic flow. Let {fj^,} denote 
a full density ergodic sequence. Then for all e < cq, 

-^[Z^c ] -ddy/p, weakly m V'^^'^\BIM). 

The proof consists of three ingredients: 

(1) By the Poincare-Lelong formula, [Z^c] = idd\og\ip'^\. This reduces the theorem to 
determining the hmit of j log 

(2) ^log|v9^| is a sequence of PSH functions which are uniformly bounded above by 
y/p. By a standard compactness theorem, the sequence is pre-compact in L^: every 
sequence from the family has an convergent subsequence. 

(3) Iv^^P, when properly normalized on each dMj. is a quantum ergodic sequence on 
OMt-. This property implies that the norm of Iv^^P on dQ is asymtotically y/p. 

(4) Ergodicity and the calculation of the norm imply that the only possible limit 
of J log Iv^l^l. This concludes the proof. 

We note that the first two steps are valid on any real analytic {M,g). The difference is 
that the norms of may depend on the subsequence and can often not equal y/p. That 
is, ^Iv^aI behaves like the maximal PSH function in the ergodic case, but not in general. 
For instance, on a flat torus, the complex zero sets of ladders of eigenfunctions concentrate 
on a real hypersurface in Mq. This may be seen from the complexified real eigenfunctions 
sin(/c,x + iC,), which vanish if and only if {k,x) G 27rZ and (k,^) = 0. Here, k G is a 
lattice point. The exact limit distribution depends on which ray or ladder of lattice points 
one takes in the limit. The result reflects the quantum integrability of the flat torus, and a 
similar (but more complicated) description of the zeros exists in all quantum integrable cases. 
The fact that ^log \ ip'^\ is pre-compact on a Grauert tube of any real analytic Riemannian 
manifold confirms the upper bound on complex nodal hypersurface volumes. 

10. Intersections of nodal sets and gedoesics on real analytic surfaces 

In ^ we discussed upper bounds on the number of intersection points of the nodal set 
with the bounary of a real analytic plane domain and more general 'good' analytic curves. 
In this section, we discuss work in progress on intersections of nodal sets and geodesies on 
surfaces with ergodic geodesic flow. Of course, the results are only tentative but it seems 
worthwhile at this point in time to explain the role of ergodicity in obtaining lower bounds 
and asymptotics. We restrict to geodesic curves because they have rather special properties 
that makes the analysis somewhat different than for more general curves such as distance 
circles. The dimensional restriction is due to the fact that the results are partly based 
on the quantum ergodic restriction theorems of |TZ2l ITZ3j . which concern restrictions of 
eigenfunctions to hyp ersurf aces. Nodal sets and geodesies have complementary dimensions 
and intersect in points, and therefore it makes sense to count the number of intersections. 

We fix {x,0 G S*M and let 

(86) 7,,5:M^Af, 7^,^(0) = x, j',^^{0) = ^ e T^M 

denote the corresponding parametrized geodesic. Our goal is to determine the asymptotic 
distribution of intersection points of jx,^ with the nodal set of a highly eigenfunction. As 
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usual, we cannot cope with this problem in the real domain and therefore analytically con- 
tinue it to the complex domain. Thus, we consider the intersections 

of the complex nodal set with the (image of the) complexification of a generic geodesic If 

(87) S, = {{t + iT eC:\T\<e} 
then 7^.^^ admits an analytic continuation 

(88) 7^,« : 5. M,. 
In other words, we consider the zeros of the pullback, 

{7:,^^^ = 0} C 5,. 
We encode the discrete set by the measure 

(89) Wf;']= Yl 

We would like to show that for generic geodesies, the complex zeros on the complexified 
geodesic condense on the real points and become uniformly distributed with respect to arc- 
length. This does not always occur: as in our discussion of QER theorems, if 7^^ ,t is the 
fixed point set of an isometric involution, then "odd" eigenfunctions under the involution 
will vanish on the geodesic. The additional hypothesis is that QER holds for 73^.^^, i.e. that 



Theorem 5.6 is valid. The following conjecture appears to be proved (|Z3j). but to be 



conservative, we state it here only conjecture: 

Conjecture 10.1. Let {M'^,g) be a real analytic Riemannian surface with ergodic geodesic 
flow. Let satisfy the QER hypothesis. Then there exists a subsequence of eigenvalues Xj^ 
of density one such that for any f G Cc{S^), 



lim Y + = / /(^)^^- 



In other words, 



k 



weak* lim ^[A^^*] = 6^-. 



in the sense of weak* convergence on Cc{S^). Thus, the complex nodal set intersects the 
(parametrized) complexified geodesic in a discrete set which is asymptotically (as A — )■ 00) 
concentrated along the real geodesic with respect to its arclength. 

This concentration- equidistribution result is a 'restricted' version of the result of ^ As 
noted there, the limit distribution of complex nodal sets in the ergodic case is a singular 
current dd^^. The motivation for restricting to geodesies is that restriction magnifies the 
singularity of this current. In the case of a geodesic, the singularity is magnified to a delta- 
function; for other curves there is additionally a smooth background measure. 

The assumption of ergodicity is crucial. For instance, in the case of a flat torus, say M^/L 
where L C is a generic lattice, the real eigenfunctions are cos(A, x), sin(A, x) where A G L*, 
the dual lattice, with eigenvalue — |Ap. Consider a geodesic 7x,^(^) = a; -|- t^. Due to the 



42 



STEVE ZELDITCH 



flatness, the restriction sin(A, xq + t^o) of the eigenfunction to a geodesic is an eigenfunction 
of the Laplacian — ^ of submanifold metric along the geodesic with eigenvalue — (A,^o)^- 
The complexification of the restricted eigenfunction is sin(A, Xq + + «t)^o) I cind its exponent 
of its growth is t|(|^,^o)|, which can have a wide range of values as the eigenvalue moves 

along different rays in L*. The limit current is idd applied to the limit and thus also has 
many limits 

The proof involves several new principles which played no role in the global result of ^ 
and which are specific to geodesies. However, the first steps in the proof are the same as in 
the global case. By the Poincare-Lelong formula, we may express the current of summation 
over the intersection points in (89) in the form, 

(90) 



[A/"; 



iddt+ir log 



Thus, the main point of the proof is to determine the asymptotics of y- log 
When we freeze r we put 

(91) ll,i{t)=^l^{t + tT). 

Proposition 10.2. (Growth saturation) If {(fj^} satisfies QER along any arcs of'^^,^, then 
in Ll^^{Sr), we have 

1 2 
lim — log 11}^%^ (t + tr) =\r\. 



Proposition 10.2 immediately implies Theorem 10.1 since we can apply dd to the 



convergent sequence log 



l*x,i'P%, (t + to obtain dd\T\. 



The upper bound in Proposition 1 10.2 1 follows immediately from the known global estimate 

1 



lim 

k—¥oo X, 



log|¥',.(7.^^(C)l<|r| 



on all of dM^-. Hence the difficult point is to prove that this growth rate is actually obtained 
upon restriction to 7^^. This requires new kinds of arguments related to the QER theorem. 

• Complexifications of restrictions of eigenfunctions to geodesies have incommensurate 
Fourier modes, i.e. higher modes are exponentially larger than lower modes. 

• The quantum ergodic restriction theorem in the real domain shows that the Fourier 
coefficients of the top allowed modes are 'large' (i.e. as large as the lower modes). 
Consequently, the norms of the complexified eigenfunctions along arcs of 7^^ 



achieve the lower bound of Proposition 10.2 



• Invariance of Wigner measures along the geodesic flow implies that the Wigner mea- 
sures of restrictions of complexified eigenfunctions to complexified geodesies should 
tend to constant multiples of Lebesgue measures dt for each r > 0. Hence the eigen- 
functions everywhere on 7^^ achieve the growth rate of the norms. 

These principles are most easily understood in the case of periodic geodesies. We let 
7x,g : 5*^ — 7- M parametrize the geodesic with arc-length (where = M/LZ where L is the 
length of 7^,5). 

First, we use Theorem |5.6| to prove 
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Lemma 10.3. Assume that {(fj} satsifies QER along the periodic geodesic 'jx,^- Let \ |7"^¥'j'| Ii2(5i-) 
be the L'^-norm of the complexified restriction of (pj along 7^^. Then, 



lim — 

A,— >oo A, 



bl I -r* (T I 1 2 II 



To prove Lemma 10.3 we study the orbital Fourier series of 7"^V'j and of its complexifi- 
cation. The orbital Fourier coefficients are 



1 

T 



VXjilx,dt))e dt, 



7 -'0 



and the orbital Fourier series is 

(92) ¥'a,(7x-,€W) = E'^a/H^ 



ne2 



Hence the analytic continuation of is given by 



(93) 



^1 



(7x,c(t + ^r)) = ^ 



ulf{n)e 



27r in (t + ir) 



By the Paley- Wiener theorem for Fourier series, the series converges absolutely and uniformly 
for |r| < eo- By "energy localization" only the modes with \n\ < Xj contribute substantially 
to the norm. We then observe that the Fourier modes decouple, since they have different 
exponential growth rates. We use the QER hypothesis in the following way: 

Lemma 10.4. Suppose that {fx^} is QER along the periodic geodesic '-)x,^- Then for all 
e > 0, there exists > so that 



E 

n:\n\>(l-e)\j 

Lemma 10.4 implies Lemma [10.3 since it implies that for any e > 0, 

n:\n\>{l-t)Xj 



To go from asymptotics of norms of restrictions to Proposition 10.2 we then use the 
third principle: 

Proposition 10.5. (Lebesgue limits) If Yx,i'^j 7^ (identically), then for all t > the 

sequence 



is QUE with limit measure given by normalized Lebesgue measure on . 



The proof of Proposition 10.2 is completed by combining Lemma 10.3 and Proposition 



10.5 Conjecture 10.1 follows easily from Proposition 6.1[ 



The proof for non-periodic geodesies is considerably more involved, since one cannot use 
Fourier analysis in quite the same way. 
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11. Nodal and critical sets of Riemannian random waves 

We mentioned above that Riemannian random waves provide a probabilistic model that 
is conjectured to predict the behavior of eigenfunctions when the geodesic flow of (M, g) is 
ergodic. In this section, we define the model precisely as in |Z4j (see also |Nic] for a similar 
model) and survey some of the current results and conjectures. We should emphasize that 
some of the rigorous results on zeros or critical points of Riemannian random waves, both 
in the real and complex domain, are much simpler than for individual eigenfuntions, and 
therefore do not provide much guidance on how to prove results for an orthonormal basis 
of eigenfunctions. But the relative simplicity of random waves and their value as predictors 
provide the motivation for studying random waves. And there are many hopelessly difficult 
problems on random waves as well, which we will survey in this section. 

For expository simplicity we assume that the geodesic flow G* of (M, g) is of one of the 
following two types: 

(1) aperiodic: The Liouville measure of the closed orbits of G*, i.e. the set of vectors 
lying on closed geodesies, is zero; or 

(2) periodic = Zoll: = id for some T > 0; henceforth T denotes the minimal period. 
The common Morse index of the T-periodic geodesies will be denoted by /3. 

In the real analytic case, {M,g) is automatically one of these two types, since a positive 
measure of closed geodesies implies that all geodesies are closed. The two-term Weyl laws 
counting eigenvalues of y/ A are very different in these two cases. 

(1) In the aperiodic case, Ivrii's two term Weyl law states 

N{X) = #{j : A, <\} = Cm Vol{M,g) + o(A'"-i) 

where m = dimM and where Cm is a universal constant. 

(2) In the periodic case, the spectrum of a/A is a union of eigenvalue clusters Cat of the 
form 

Cn = {(^)(A^ + + ^^'^ i = l...dN} 

with fim = 0{N^^). The number c^at of eigenvalues in Ca? is a polynomial of degree 
m — 1. 

We refer to |HoI-IVt IZ4j for background and further discussion. 

To define Riemannian random waves, we partition the spectrum of a/A^ into certain 
intervals In of width one and denote by Hjj^ the spectral projections for a/"A^ corresponding 
to the interval In. The choice of the intervals In is rather arbitrary for aperiodic {M,g) and 
as mentioned above we assume In = [N,N + 1]. In the Zoll case, we center the intervals 
around the center points '^N + ^ of the A^th cluster Cn- We call call such a choice of 
intervals a cluster decomposition. We denote hj dN the number of eigenvalues in In and put 
Hn = ranll/^ (the range of H/^^). 

We choose an orthonormal basis {^Nj}'j=i for T-Ln- For instance, on S'^ one can choose the 
real and imaginary parts of the standard Y^'s. We endow the real vector space T-Ln with 
the Gaussian probability measure 7Ar defined by 

(7 \ djv/2 dx 
e-"-l^l'dc, f = J2cjipNj, dN = dimUN. 
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Here, dc is ciAr-dimensional real Lebesgue measure. The normalization is chosen so that 
E^jy (/, /) = 1, where E^^^ is the expected value with respect to 7Ar. Equivalently, the 
real variables Cj {j = l,...,d]y) are independent identically distributed (i.i.d.) random 
variables with mean and variance 7^5—; i.e., 

We note that the Gaussian ensemble is equivalent to picking /jv G "Hjv at random from the 
unit sphere inUN with respect to the inner product. 

Depending on the choice of intervals, we obtain the following special ensembles: 

• The asymptotically fixed frequency ensemble T-Li^, where Ix — [A, A + 1] and where 
T-Li^ is the vector space of linear combinations 

(95) /a= Y1 

j:A,e[A,A+l] 

of eigenfunctions with A^ (the frequency) in an interval [A, A + 1] of fixed width. 
(Note that it is the square root of the eigenvalue of A, not the eigenvalue, which is 
asymptotically fixed). 

• The high frequency cut- ojf ensembles ^[o,a] where the frequency is cut-off at A: 

(96) /a= E 

• The cut-off Gaussian free field, 

(97) /a = E ^= f"- 

i:A,<A ^ 

One could use more general weights w{\j) on a Sobolev space of functions or distributions 
on M. In the physics terminology, w{\j) (or it square) is referred to as the power spectrum. 

The key reason why we can study the limit distribution of nodal sets in this ensemble is 
that the covariance kernel 



(98) n,,(a;,y)=E^,(/jv(x)/^(y))= ^ Vx,{x)vxM^ 

j:\jelN 

is the spectral projections kernel for y/A. 

11.1. Equidistribution of nodal sets for almost all sequences of random waves. 

The real zeros are straightforward to define. For each fx G Ti[o,x] or Hi^ we associated to 
the zero set Zf^ — {x e M : fx{x) — 0} the positive measure 

(99) {\Zf,\,i^)= [ ^dn^-\ 

where dl-V^'^ is the induced (Hausdorff) hypersurface measure. 
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The main result we review is the hmit law for random sequences of random real Riemannian 
waves. By a random sequence, we mean an element of the product probability space 

oo N 

(100) -Hoo = n = n 

N=l N=l 

Theorem 11.1. |Z4j Let {M,g) be a compact Riemannian manifold, and let {/^} be a 



random sequence in (100). Then 



1 1 



— > — |%J ^ dVg almost surely w.r.t. {H oo 1 Too ) ■ 



n=l 

11.2. Mean and variance. We first show that the normalized expected limit distribution 
l-E|Zj^| of zeros of random Riemannian waves tends to the volume form dVg as A — )■ oo. 
That is, we define the 'linear statistic', 

(101) X^r(/;v) = (^,|^/,|), ^gC(M) 
and then define 

(102) (E,^|Zy^U)=E,^Xir, 

Theorem 11.2. Let {M,g) be a compact Riemannian manifold, let Hio^x] be the cutoff ensem- 
ble and let {T-LniIn) be the ensemble of Riemannian waves of asymptotically fixed frequency. 
Then in either ensemble: 

(1) For any {M,g), lim^^^ j^E^^{\ZfJ,^) = j^^^dVg. 

(2) For a real analytic {M,g), Var{^X^)) < C. 

We restrict to real analytic metrics in (2) for the sake of brevity. In that case, the variance 



estimate follows easily from Theorem 2.1 



11.3. Density of real zeros. The formula for the density of zeros of random elements of 
% N can be derived from the general Kac-Rice formula |BSZll IBSZ2t INicj : 

(103) nZf^\=K^{z)dVg, K^ix) = J DiO,^,x)m\d^. 

Here, D[q,^,x)dqdC, is the joint probability distribution of the Gaussian random variables 
(■^(x), V'0(a;)), i.e. the pushforward of the Gaussian measure on Tix) under the map — )■ 
lijlx),Vijlx)). Note that the factor det(^D in |B'Hzn[BgZ2i equals in the codimension 
one case. Indeed, let df* be the adjoint map with respect to the inner product g on T^M. 
Let dfx o df* : M — i- M be the composition. By det dfx o df* is meant the determinant with 
respect to the inner product on T^M; it clearly equals in the codimension one case. 
The formulae of jBSZlt IBSZ2j (the 'Kac-Rice' formulae) give that 

(104) D{0,^-z) = Z^{z)D^{^;z), 
where 

(105) D^iC,z) = ^^exp(-(A-ie,0) 

7r™vdet A 



EIGENFUNCTIONS AND NODAL SETS 47 

is the Gaussian density with covariance matrix 

(106) K = C- B*A-^B = {CI, - BqA-^Bq) , (g = 1, . . . , m) 
and 

(107) Z(x- ^ ' 



VdetA ny/A ' 



Here, 



A = A^ix 



q,q' = 1,. . . ,m. 



Making a simple change of variables in the integral (103), we have 

Proposition 11.3. |BSZlj On a real Riemannian manifold of dimension m, the density of 
zeros of a random Riemannian wave is 

(108) ^^■''"(^) = .~lV.,'n,„,„, -<»■■■ ll^"'^)'^'^ll^-P<-«-^»'^- 

where A^(x) is a symmetric form on T^M. For the asymptotically fixed freqency ensembles, 
it is given by 

If 1 

^^(^) = zr {^^'^ dyUij^{x,y)\^=y - - — -d^Ili^{x,y)\^=y (g) dyUi^{x,y)l 



In the cutoff ensemble the formula is the same except that Uj^^ is replaced by II[o,n] ■ 

We then need the asymptotics of the matrix elements of A^(a;). They are simplest for the 
round sphere, so we state them first in that case: 

Proposition 11.4. LetllN : L'^{S'^) Hn be the orthogonal projection. Then: 

• (A) Un{x,x) = voiis"-) ^^' 

• (B) d^UN{x,y)\^=y = dynN{x,y)\^=y = 0; 

• (C) d^® dyIlN{x,y)\x=y = mVol{S"') ^^NdN9x- 

We refer to |Z4j for the calculation, which is quite simple because of the invariance under 
rotations. The expected density of random nodal hypersurfaces is given as follows 

Proposition 11.5. In the case of S"^, 

(109) K^{x) = CmXN-^C^N, 
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where Cm = ^ /r™ |^| exp (-(^, ^)) d^. 



Proof. By Propositiosn 11.3 and 11.4, we have 

^/Vol{S^ 



(110) 
where 



TX' 



dN \mVol{S''' 



A^(x)m|exp(-(e,0)rfe, 



-X^d^gx 



□ 



11.4. Random Riemannian waves: proof of Theorem 11.2 We now generahze the 
result to any compact C°° Riemannian manifold (M, g) which is either aperiodic or Zoll. 
As in the case of 5*™, the key issue is the asymptotic behavior of derivatives of the spectral 
projections 

(111) n7^(x,?/)= ^ (^A,(a;)^A,(y). 

Proposition 11.6. Assume {M,g) is either aperiodic and In = [A^, + 1] or Zoll and In 
is a cluster decomposition. Let II/^ : L'^{M) Hn be the orthogonal projection. Then: 

• (A) n7^(x,x) = ^^j/^^^^ dN{l + o{l)); 

• (B) 4n/^(a;,: 



Vol{M,g))'^N^ 

\x=y = dyUN{x,y)\x=y 



• (C) d^®dylii^{x,y%=y = y^^^j^^XjjdNgxil + o{l)). 
In the aperiodic case, 
(1) n[o,A](x,x) = C„A- + o(A— 1); 



(2) 4. o rfj^n[o,A]i 



X 



, y )\x=y 



CmX^+^gx + o(A 



m+l^ 



In the Zoll case, one adds the complete asymptotic expansions for 11/^ over the N clusters 
to obtain expansions for Hn- 

We then have: 

Proposition 11.7. For the asymptotically fixed frequency ensemble, and for any {M,g) 



which is either Zoll or aperiodic (and with In as in Proposition 11.6) , we have 



(112) 



K^x) = ;^^^.(X^/M™lie||exp(-^(e,e))^ 



where Cm = ^ /^m ||^|| exp {—{C,,^)) d^. The same formula holds for the cutoff ensemble. 

Proof. Both on a sphere S"^ o r on a mo re general {M,g) which is either Z oll or aperiodic, 
we have by Propositions 11.4 resp. 11.6 and the general formula for in [ 11.3 that 



Vol{M,g) 



1 + 0(1)) 0(1) 

0(1) N'g^{l + o{l)) 



(113) 
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It follows that 

(114) = - B^*{A^)-'B'' = y^j^f' 9. + o{N). 

Thus, we have 

K^^x) ~ ^^^^P^/«J|A^(x)n||exp(-(e,0)c^e 

(115) 

= ^ A. ||(/ + o(l))(x)m| exp (-(e, 0) rfe, 

where o(l) denotes a matrix whose norm is o(l), as as — i- oo we obtain the stated 
asymptotics. 

□ 

So far, we have only determined the expected values of the nodal hypersurface measures. 
To complete the proof of Theorem 11. 2[ we need to prove: 



Proposition 11.8. If {M,g) is real analytic, then the variance of -^X^ is hounded. 



for /iv G "H/jy, j^Zfj^ has bounded mass. Hence, the random 

□ 



Proof. By Theorem 2.1 

variable -^X^ is bounded, and therefore so is its variance. 



Remark: 

The variance of -^X^ is given by 

(116) Var{^X^) = ^[ [ {K^ {x,y) - {x)K^ {y)) ^lj{x)^l,{y)dVg{x)dV,{y), 

where i^^(x,?/) = E^j^(Z/^(x) ® Zfj^iy)) is the pair correlation function for zeros. Hence, 
boundedness would follow from 

(117) / [ K^{x,y)dV,{x)dVg{y)<C. 

J M J M 



There is a formula similar to that for the density in Proposition 11.3 for K2{x,y) and it is 
likely that it could be used to prove boundedness of the variance for any C°° Riemannian 
manifold. 



11.5. Random sequences and proof of Theorem 11. 1[ We recall that the set of ran- 



dom sequences of Riemannian waves of increasing frequency is the probability space Tioo = 
Y[n=i^In with the measure 7oo = n?^=i7A'- element in T-Loo will be denoted f = {/at}. 
We have. 

An an 

By a density argument it suffices to prove that the linear statistics Y^{Zfj^,ilj) — y^^j^^ fj^ ipdVg — )■ 
almost surely in Tioo We know that 

(i) lim^^^ ^ Ek<N^it^^) = y^^lfe) iM'I'dVg; 

(ii) Var^Y^X^) is bounded on Hoo- 
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Since -^X^ for {, = 1, 2, . . .} is a sequence of independent random variables in "He 
with bounded variances, the Kolmogorov strong law of large numbers implies that 



almost surely. 



11.6. Complex zeros of random waves. We now state a complex analogue of the equidis- 
tribution of real nodal sets and show that it agrees with the the limit formula of Theorem 

EH 

We complexify Riemannian random waves as 

rC \ ^ C 

Jn — / ^'^Nj^Nj- 
J=l 

We note that the coefficients cjsij are real and that the Gaussian measure on the coefficients 
remains the real Gaussian measure 7iv- The two point function is the analytic extensions to 
the totally real anti-diagonal in Mq x Mq is therefore 

(118) E(|/^(C)p) = n,,(C,C)= E l^5'(0r- 



As in the proof of Theorem 9.1, the current of integration over the complex zero set 

Zfc = {C e Mc : = 0} 

is the (1, 1) current defined by 

' N 

for smooth test forms of bi-degree (m — l,m — 1). In terms of scalar functions if) we may 
define ZfC as the measure, 

m— 1 



where Ug = iddp is the Kahlermetric adapted to g. 

The proof of the next result is close to the proof of Theorem and we therefore refer to 
[Z4j for the details: 

Theorem 11.9. |Z4] Let {M,g) be a real analytic compact Riemannian manifold. Then for 
either of the Riemannian random wave ensembles 

E,, (^[%]) ^ ^dd\^\,, weakly zn V'^-^'^^B^^M). 

As mentioned above, this result shows that the complex zeros of the random waves have 
the same expected limit distribution found in |Z3j for real analytic compact Riemannian 
manifolds with ergodic geodesic flow. 
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12. Percolation heuristics 

In this final section, we review some of tfie more speculative conjectures relating nodal 
sets of both eigenfunctions and random waves to percolation theory. The conjectures are 
often quoted and it therefore seems worthwhile to try to state them precisely. The only 
rigorous result to date regarding eigenfunctions is the theorem of Nazarov-Sodin on the 
expected number of nodal domains for random spherical harmonics |NSj (see jZ5j for a brief 
over- view) . 

The percolation conjectures concern the statistics of sizes of nodal domains or nodal com- 
ponents. They are based on the idea that the nodal domains resemble percolation clusters. 
One might measure the 'size' of a nodal component Ax. by its hypersurface area 'H'^~^{Axj), 
and a nodal domain Dxj by its volume 'H"(-Daj) • Let us restrict to the case of surfaces. 
For the purposes of this article, we introduce the term length spectrum of the nodal set as 
the set 

(119) Lspi^x) = m\Cx,) : Z^, = [jCx;} 

of lengths of its components, counted with multiplicity. It is encoded by the empirical 
measure of surface areas 

(120) rf/iL = E'^^HC.,,) e Pi(M), 

(where V{fl) is the set of probability measures on Q), or equivalently by the length distri- 
bution function, 

(121) Cx{t)= E ^'(Ca,). 
We also consider the area spectrum, 

(122) Aspicpx) = m\Ax,) : M\Z^^ = [JAx.,}, 

encoded by its empirical measure It is encoded by the empirical measure of surface areas 

(123) *.^ = a;4?)^''"'-'«'"^''*>' 

or by the area distribution function, 

(124) Ax{t)= E ^'(^A.^)- 

Of course, there are some obvious constraints on such spectra; e.g. in the analytic case, there 
could only exist 0(A) components with "H^-length of order 1, and only a bounded number 
of order A. 

In computer graphics of eigenfunctions on plane domains or surfaces, one sees many 'small' 
components Cxj of the nodal set whose length appears to be of order j. But one also sees long 
snaky nodal lines. How long are they? Do they persist as A — )■ oo? Roughly speaking, one 
may ask what proportion of the components come in sizes with different orders of magnitude. 
Of course, this depends on how many components there are, so it could be simpler to work 
with C{Lpx),A{ipx)- 
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• How many components have "H"" ^-surface measure which is > CX* for some given 
< 7 < 1. It is possible that some individual nodal component has ^/""^-surface area 
commensurate with that of the entire nodal set, as in the Lewy spherical harmonics 
with just two or three nodal components [Lew]. 

• How many components have ^{""^-surface measure (i.e. length in dimension two) 
which is bounded below by a constant C > independent of A? Such components 
are sometimes termed "percolating nodal lines" since their hypersurface volume is 
commensurate with the size of the macroscopic object (i.e. M). 

• How many components have ?/'^~^-surface measure of the minimal order ^? 

The percolation conjectures relate the asymptotic distribution of lengths of nodal com- 
ponents and areas of nodal domains of eigenfunctions as defined in (120)-(123) to lengths 
of boundaries and areas of percolation clusters at criticality. There are different types of 



conjectures for the fixed frequency ensemble and the high frequency cutoj^ ensemble (see { 11 



for the definitions). According to the random wave hypothesis, the conjectures concerning 
the fixed frequency ensemble (e.g. random spherical harmonics of fixed degree) should also 
apply to nodal sets of eigenfunctions of quantum chaotic systems. 

Percolation theory is concerned with connectivity and transport in a complex system. In 
particular, it studies connected clusters of objects in a random graph. In bond percolation 
the edges of the graph are independently open or closed with some probability p. The open 
edges form a subgraph whose connected components form the clusters. In site percolation 
the vertices are open or closed and an open path is a path through open vertices. The open 
cluster C{v) of a vertex is the set of all open vertices which are connected to v by an open 
path. 

There also exists an analogous continuum percolation theory for level sets of random 
functions. We will assume the random functions are Gaussian Riemannian random waves 
on a surface. The main problem is to study the connectivity properties of level sets {/ = t}. 
One imagines a random landscape of lakes and islands depending on the variable height t 
of the water, the islands being the super-level sets {/ > t} of the random functions. For 
high water levels, the islands are disconnected, but as the water level is lowered the islands 
become more connected. At a critical level they 'percolate', i.e. it is possible to traverse 
the landscape while remaining on the land. A review with many illustrations is given by 
Isichenko [lii] (see Section E (c), pages 980-984). As explained in [lii] page 984, the contour 
lines of a random potential are associated to hulls of percolation clusters. Hence the area 
spectrum (122) is similar to the set of sizes of connected clusters in a percolation model. 

In the physics literature, the random functions are usually functions on M? (or possibly 
higher dimensional M") and the Gaussian measure on the space of functions corresponds to 
a Hilbert space inner product. The Hilbert space is usually taken to be a Sobolev space, 
so that the inner product has the form J w{S,)\f{C)\^d^ (where / is the Fourier transform 
of /) and w(^) = The case C = is known as the Gaussian free field (or massless 

scalar field) and is quite special in two dimensions since then the inner product /j^a |V/p(ix 
is conformally invariant. There are rigorous results on level sets of discretizations of the 
Gaussian free field and their continuum limits in |SSl IMij . with authoritative comments on 
the physics literature. 

For purposes of this exposition, we assume the Riemannian random waves fall are of the 



types discussed in {11 In all cases, we truncate the frequency above a spectral parameter A 
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and consider asymptotics as X oo. In this high frequency hmit, the random waves oscillate 
more rapidly on the length scale j. Since the conjectures and results depend strongly 
on the chosen weight w, we break up the discussion into two cases as in ^ pTj the high 
frequency cutoff ensemble and the fixed frequency ensemble. For each ensemble we let Ea 
denote the expectation with respect to the Gaussian measure on the relevant space of linear 
combinations. Then we may ask for the asymptotic behavior of the expected distribution of 
lengths of nodal lines, resp. area of nodal domains 

(125) ExdfiL, ExdfiA, 



where dfii, resp. dfiA are the empirical measures of lengths (120) of nodal lines, resp. areas 



(123) of nodal domains. 



12.1. High frequency cutoff ensembles. The distribution of contour lengths of certain 
Gaussian random surfaces over was studied at the physics level of rigor in jKH] . They 
define the Gaussian measure as e'^^^^^^dh where the 'free energy' is defined by 

where x is a cutoff function to [0, 1] (they use the notation a for \ in our notation). When 
C = 0, this is a truncated Gaussian free field (truncated at frequencies < A) and its analogue 
on a surface (M, g) is the Riemannian random wave model with spectral interval [0, A] and 
weight w{X) = J. The parameter ( is referred to as the 'roughness exponent' in the physics 
literature. In the case of the Gaussian free field C = the inner product is the Dirichlet 
inner product J^2 |V/p(ix. 

An important feature of the ensembles is scale-invar iance. In the special case ( = 
(and dimension two), the Dirichlet inner product |V/|g(iylg is conformally invariant, i.e. 
invariant under conformal changes g — )■ e'^g of the Riemannian metric. When ( ^ this is 
not the case, but it is assumed in |KHj that the fluctuations of the random Gaussian surface 
with height function h are invariant under the rescaling h{r) — > c^'^h{cr) for any c > 1. The 
authors of |KHj then make a number of conjectures concerning the distribution of contour 
lengths, which we interpret as conjectures concerning Kdfii- First, they consider contours 
(i.e. level sets) through a fixed point xq and measure its length with the re-scaled arc-length 
measure Xds, i.e. with arclength s in units of j. They define the fractal dimension of a 
nodal line component as the dimension D so that s ~ where R is the radius of the 
nodal component (i.e. half the diameter). They define P{s) as the probability density that 
the contour through xq has length s. The principal claim is that P{s) ~ satisfies 
a power law for some exponent r ( |KHj (4)). They also defines the distribution of loop 
lengths P{s) ~ P{s)/s as the probability density that a random component has length s. 
We interpret their P{s) as the density of hm^^-oo Ed/i^ with respect to ds on M. We thus 
interpret their conjecture as saying that a unique weak* limit of this family of measures 
exists and has a density relative to ds with a power law decay as above. 

The claims are based in part on scaling properties of the contour ensemble. They also 
are based in part on the expectation that, at 'criticality', the key percolation 'exponents' of 
power laws are universal and therefore should be the same for the discrete and continuum 
percolation theories (see e.g. |IsiK] ). In |KHj . the authors suggest that when a certain 
roughness exponent ( vanishes (the critical models), the continuum problem is related to 
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the four-state Potts model. The g-state Potts model is an Ising type spin model on a lattice 
where each spin can take one of q values. It is known to be related to connectivity and 
percolation problems on a graph |Baxl IWuj . 

They compute D,t by relating both to another exponent xi defined by a "contour cor- 
relation function" Qi{r), which measures the probability that points at x,x + r lie on the 
same contour loop. They claim that Qi{r) ~ |r|~^^^ They claim that D(3 — r) = 2 — 2xi 
and D{t — 1) = 2 — ^. As a result, D = 2 — xi — C,/2,t — 1 = 2~-l^-c/2 ■ From the mapping 
to the four-state Potts model, they conclude that xi = |. 

There exist rigorous results in (SS[ IMij relating discretizations of the Gaussian free field 
(rather than high frequency truncations) to the percolation models. They prove that in 
various senses, the zero set of the discrete Gaussian free field tends to an SLE^ curve. It 
does not seem to be known at present if zero sets of the high frequency truncation of the 
Gaussian free field also tends in the same sense to an SLE4 curve. Note that the SLE curves 
are interfaces and that one must select one component of the zero set that should tend to 
an SLE curve. There might exist modified conjectures regarding CLE curves. 

To determine the 'critical exponents' in continuum percolation, it is tempting to find a 
way to 'map' the continuum problem to a discrete percolation model. A geometric 'map' 
from a random wave to a graph is to associate to the random function its Morse-Smale 



decomposition, known in the physics literature as the "Morse skeleton" (see ^2.6 or |Web] 



for an extensive exposition). As discussed in |Weij . and as illustrated in Figure 10 of [Hi] 
the Morse complex of the random function plays the role of the lattice in lattice percolation 
theory. 



12.2. Fixed frequency ensembles. We now consider Riemannian random waves of asym- 
totically fixed frequency A, such as random spherical harmonics of fixed degree or Euclidean 
random plane waves of fixed eigenvalue. In this case the weight is a delta function at the 
frequency. One would expect different behavior in the level sets since only one frequency is 
involved rather than the superposition of waves of all frequencies < A. 

A recent exposition in the specific setting of random Euclidean eigenfunctions of fixed 
frequency is given by jEGJSj . The level sets play the role of open paths. Super-level sets are 
compared to clusters of sites in a critical 2D percolation model, such as bond percolation 
on a lattice. Each site may of the percolation model may be visualized as a disc of area 
i.e. as a small component. The nodal domains may be thought of as connected clusters 
of a number n such discs. Since nodal domains are connected components in which the 
eigenfunction is either positive + or negative — , they are analogous to clusters of 'open' or 
'closed' vertices. 

The main conjectures in this fixed frequency ensemble are due to E. Bogolmony and C. 
Schmidt |BSj . They conjecture that the continuum percolation problem should belong to 
the same universality class as the Potts model at a certain critical point (where q is related 
to a certain temperature) for a large rectangular lattice and that the nodal lines in the 
A — )■ cxD limit tend to SLEq curves. This is similar to the predictions of |KH] but for a very 
different ensemble where there is little apriori reason to expect conformal invariance in the 
limit. There are parallel conjectures in |BBCF] for zero-vorticity isolines in 2D turbulence, 
which are also conjectured to tend to SLEq curves. They remark (page 127) that this limit 
is surprising since continuous percolation models assume short-correlations in the height 
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functions whereas the vorticity field correlations decay only like r""^/^. They write, "When 
the pair correlation function falls off slower than r~^/^, the system is not expected generally to 
belong to the universality class of uncorrelated percolation and to be conformally invariant" . 
The same remarks apply to the fixed frequency ensemble, where the correlation function is 
the spectral projection n[A,A+i](a^, y) for a fixed frequency. In this case, the correlations decay 
quite slowly as r~5; we refer to |BS2j for this background and also for an argument why the 
nodal sets should nevertheless resemble conformally invariant SLEg curves. 

If the nodal lines in the fixed frequency model are equivalent to the critical percolation 
model, then the 'probability' of finding a nodal domain of area s should decay like s~'^ where 
r = ^ > 2 (see |SAj . p. 52 for the percolation theory result). Under some shape assump- 
tions adopted in |EGJSj . it is equivalent that the probabihty of finding clusters consisting 
of n discs is of order n~'^. For random spherical harmonics, one may ask for the probability 
that a spherical harmonic of degree N has size n. For a fixed (M, g) with simple eigenvalues, 
this notion of probability from percolation theory does not make sense, but we might assume 
that the number of of nodal components is of order and ask what proportion of the nodal 
components has size 1. To obtain a percolating nodal line, one would need a cluster with 
n = X sites, and thus the proportion of such nodal components to the total number would 
be of order A""^. Thus, if there are CX^ total components, the number of such components 
would be around A^""^ = X~^ < 1, so the model seems to predict that such macroscopic 
nodal lines are quite rare. It also predicts that the 'vast majority' of nodal components are 
close to the minimal size, which does not seem so evident from the computer graphics. 
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